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> 

k><( I The application of the effective field theory (EFT) method to nuclear systems 

^H ' is reviewed. The roles of degrees of freedom, QCD symmetries, power count- 

ing, renormalization, and potentials are discussed. EFTs are constructed for 
the various energy regimes of relevance in nuclear physics, and are used in 
systematic expansions to derive nuclear forces in terms of a number of pa- 
rameters that embody information about QCD dynamics. Two-, three-, and 
many-nucleon systems, including external probes, are considered. 
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I. INTRODUCTION 

There is little doubt that QCD is the correct theory of the strong forces that bind quarks 
into baryons and mesons and yet, despite many efforts, it has not been possible to derive the 
interaction among these composite states from first principles. At the root of the problem 
is our inability to solve the dynamics of QCD at hadronic scales because of a large coupling 
constant. Much of the description of hadronic properties and interactions has been limited 
to phenomenological models whose connections with QCD are unclear. 

The strong-coupled dynamics among baryons is particularly rich: it gives rises to bound 
states, of which ordinary nuclei are particular interesting. A superficial look at hadronic and 
nuclear tables conveys complexity, but a closer look reveals some regularities that can be 
translated into different mass scales. At least two scales influence the bulk of the dynamics 
in a nucleus. The masses of the main constituents (nucleons) and their excitations (chiefly 
the delta isobar), the masses of most mesons that can be exchanged among them (the rho, 
omega, etc.), and the chiral symmetry breaking scale; they all cluster near what we can call 
the characteristic QCD scale, Mqcd ~ 1 GeV. On the other hand, the reciprocal of the 
sizes of light nuclei, the Fermi momentum of equilibrium nuclear matter, the pion mass and 
decay constant, and the delta- nucleon mass difference take values closer to a lighter scale 
Mnuc ~ 100 MeV. From these two scales, a third, MI^^^/Mqcd ~ 10 MeV, manifests itself 
in the order of magnitude of nuclear binding energies. 

In nuclear physics there is thus a natural separation of scales, and one can hope to use 
it to simplify the theoretical analysis. Such an approach has proved useful in other areas of 
physics, and is at the heart of the use of renormalization group techniques. The method of ef- 
fective field theory (EFT) evolved from these techniques, and it was designed to incorporate 
and benefit from the existence of separate scales, even when couplings are not small com- 
pared to 1. It seems to constitute an ideal tool for low-energy nuclear physics, a tool which 
has similarities but also significant differences compared to the more traditional phenomeno- 
logical approach. The method is firmly based on the known (approximate) symmetries of 
QCD, but parametrizes our ignorance of dynamical details in a number of undetermined 
constants. This article is a review of the program of research that attempts to describe 
nuclei and processes involving them using the method of EFT. 

The scales described above are consistent with a picture of nuclei as composed of non- 
relativistic nucleons of typical momenta Q ~ Mnuc and energies Q^/Mqcd- Now, slow 
particles make for poor microscopes. A nucleon in the nucleus most of the time cannot 
resolve features with a spatial extent ^ 1/Q- In particular, particles that live for short 
times and distances can in first approximation be treated as points. Only those particles 
that propagate over distances longer than ~ 1/Q need to be retained in the theory as explicit 
degrees of freedom. They will supply the correct low-energy analytic structure (poles and 
cuts) of physical amplitudes. This is a significant simplification, as it means that a large 
number of states (massive conventional hadrons or exotic quark states) can be accounted 
for just through their contributions to local interactions. Further simplification follows from 
the fact that we have been able to identify a number of symmetries of QCD, as well as the 
pattern of their breaking. Low-energy S'-matrix elements have to exhibit the same symmetry 
pattern, which markedly constrains the possible interactions between the relevant degrees 
of freedom. This still leaves us with infinitely many local interactions, as arbitrary numbers 



of derivatives appear to account for corrections to the point approximation. Worse still, 
parameters characterizing the interaction strengths need not be small. The crux of the EFT 
method is to find a way to organize this infinite set of interactions using an expansion in 
the small parameter Q. 

The method can be illustrated with a very familiar example. Let us suppose we want to 
study a system whose overall size is Ri but that contains smaller structures of size Rg <^ Ri. 
We can do so by shining light of wavelength A onto the system, but for my purposes here it 
does not really matter if this "light" comes from a laser source or from electron scattering or 
whatever; in fact I will consider just a static field, the only important thing being that there 
is a scale A associated with its (spatial) variation. We can learn about the system because 
its charge density p{r) affects the interaction energy between the system and light. 

Let us consider first the case when A ^ i?/. In classical electrodynamics we argue that 
the interaction energy can under these circumstances be written as a multipole expansion. 
In first approximation, we need only to know the overall charge q of the system and the 
electromagnetic potential at its center. In second approximation we need also the system's 
dipole moment and the gradient of the potential, and so on. Schematically, 

d^rp{f)(f){f) = g 51 C" ( -r ) • (1) 

n=0 ^ -^ ^ 

This expansion is justified by dimensional analysis: except for very pathological charge 
distributions, the interaction energy is analytic in 1/A and successive terms must come in 
powers of Ri/X times dimensionless factors c„ of 0(1). The latter encode all we need to 
know about the system; they can be calculated from p(r), but if p{r) is not known, we do 
not have to weigh on every single dimple of the charge distribution. To a given accuracy, 
we can truncate the multipole expansion and parametrize the electromagnetic properties 
of the system by a finite number of parameters instead of a continuous pir). The role 
of an exact (approximate) symmetry is that it forbids (suppresses) certain terms in the 
expansion. Now, all this can be translated to QED. All we need to do is to write the most 
general Lagrangian allowed by gauge invariance that involves fields for the system (and 
perhaps its excited states) and the photon. This Lagrangian will contain not only minimal 
couplings governed by charge but also operators with more derivatives, such as a Pauli 
term. The rules of field theory then can be used to derive the electromagnetic form factors 
of the target — or more generally any scattering amplitude — with contributions ordered by 
Ri/X. There is a difference, however. Quantum fluctuations from loops can probe small 
distances. The procedure of renormalization is designed to ensure that the physics of these 
small distances contribute to the physics of large distances only in an average sense, and 
that the multipole expansion remain valid. Left behind are some calculable cuts. Because 
(point) charge interactions are renormalizable, the non-analytic terms are suppressed only 
by powers of the fine-structure constant, a. Loops involving higher- derivative operators are 
further suppressed by powers of Ri/X. 

Consider now an increase in resolution: A ~ i?;. Again, because X ^ Rg the contributions 
of the small structures can be treated in an expansion in Rg/X as in Eq. (|l]), except that the 
multipole parameters now refer to the inner charge distribution of charge g^. Contributions 
from the outer cloud, on the other hand, can only be accounted for by the full knowledge of 
the outer distribution pi{r). The interaction energy is now 



j d\p{f)<f>{r} =q,Y: cl:^ (y ) + J d'rp,{r)<p{r). (2) 



In the quantum context these outer contributions that are non-analytic in 1/A come from 
loops generated by the light particles flying within the target at the time light strikes. The 
magnitudes of loop diagrams will depend crucially on the interactions among the structures 
that make up the system. In general, they are not small compared to tree-level graphs. 

We can translate this simple example to a typical nucleus, where Mqcd ~ 1/-Rs, 
Mnuc ~ 1/-R«; and Q ~ 1/A. It turns out that in this case we can justify not only the 
"multipole" or derivative expansion, but we can control the loop expansion as well. 

The first (and simplest) situation, Q ^ Mnuc, corresponds to nucleons of momenta suffi- 
ciently small that all interactions among nucleons can be described as contact interactions. 
One might be tempted to dismiss this as an uninteresting case, but it is not so, as some 
information on the two-nucleon system reveals. At low momenta the scattering amplitude 
in both isospin channels (/ = and / = 1) can be very well- described by the effective range 
expansion. One finds that, as expected on the basis of dimensional analysis, most of the S- 
wave effective range parameters scale with M~J^. However, the scattering lengths a2 are large 
and reflect the existence of a new, smaller mass scale which I denote by K ~ l/a2- In the ^5*1 
channel, K^^'^^^ ~ 40 MeV, somewhat smaller than Mn^c- In the ^5*0 channel, K^^'^'o) ~ 8 MeV 
^ Mnuc- For simplicity, I take an average K ~ 30 MeV. The existence of a new scale requires 
some amount of flne-tuning. whatever its origins, the consequence is that these channels 
have anomalously shallow bound states with binding energies ~ 'R'^/Mqcd ^ ^nud^QCD 
In the singlet channel this bound state is virtual, but in the triplet it is real — the deuteron. 

The idea of a low-energy expansion is in fact as old as nuclear physics: already in the 
30 's Bethe and Peierls ||l| reasoned that, because the deuteron has a large size A compared 
to the range Ri of nuclear forces, up to an error 0{Ri/X) the sole effect of the potential 
is to provide an energy-independent boundary condition at r ~ i?; to an otherwise free 
Schrodinger equation. This can be recast in terms of an EFT of non-relativistic nucleons 
— non-relativistic because Q is small compared to the nucleon mass. As we will see in 
detail below, interactions contain K besides Mnuc and the appearance of Q"^ /Mqcd as the 
only scale in energy denominators of loop diagrams generates enhancements that require a 
resummation of a large class of diagrams. After resummation, the EFT includes terms of 
all orders in Q/K but is perturbative in Q /Mnuc- Although the EFT in the two-body case 
is little more than a generalization of the technique used by Bethe and Peierls, it allows 
an easy extension of these arguments to processes where low-energy photons, electrons and 
neutrinos probe the deuteron. It also allows the study of nucleon-deuteron scattering at 
low-energies, and perhaps of the triton and other nuclei as well. I will refer to momenta 
Q < Mnuc as "very low" . 

As we increase momenta and Q ~ Mnuc-, pions can no longer be considered heavy, 
as Q is comparable to their mass. The EFT has to be augmented to include them as 
explicit degrees of freedom. (Because the delta-nucleon mass difference is of the same order 
of magnitude, explicit inclusion of the delta might also be called for.) Explicit inclusion 
of the pion dynamics gives rise to non-analytic terms analogous to those in our simple 
electrodynamics example. All short-range dynamics from massive states can be treated in 
an expansion in Q/Mqcd-, and is parametrized in analogy to the c^n^ in Eq. (0). What is 



crucial for the success of the approach is that one can show that — in a sense to be made 
more precise below — pion loops can also be ordered in terms of an expansion in Q/Mqcd- 
In terms of Eq. (^, this would mean that, apart from logarithms and dimensionless factors 
of 0(1), contributions from the outer cloud would come with a small factor {Ri/AnXy for 
each loop, namely 

Jd'rpi{r}<P{r) = q,Y.lL[-^) , (3) 

where jl are calculable functions of the Lagrangian parameters. It is chiral symmetry that 
effects this magic by requiring that pion interactions be proportional to Q, and therefore 
weak at low energies. 

The EFT appropriate to these momenta has been extensively applied to systems contain- 
ing one or no nucleon, and in that context it is called Chiral Perturbation Theory (^PT). 
The application of xPT to nuclear systems is complicated by the resummation needed to 
accommodate the infrared enhancement. What class of diagrams is to be re-summed? The 
answer to this question depends on the relative magnitude of pion and shorter-range effects. 
It is reasonable that at sufficiently small Q, pions can be treated perturbatively, while as Q 
reaches some scale Mj^^ pion effects become non-perturbative. If M^n is as low as Mnuc, 
then the theory with perturbative pions is of little use, since then the pion can be simply 
replaced by a point. If Mnn is as high as Mqcd, then non-perturbative pions are not nec- 
essary as long as Q < Mqcd f\- I will refer to momenta Q ^ M^n as "low" and to momenta 
MisfN ^ Q ^ Mqcd as "moderate". (What the scale M^n is is a matter of current debate.) 

What is unique about the application of EFT to nuclear physics is this fascinating in- 
terplay between perturbative and non-perturbative phenomena. EFTs already have a long 
history, as its first example likely is the work of Euler and Heisenberg in QED 0. However, 
its modern concept evolved only during the late 60's and the 70's, as a confluence of several 
ideas. First, in S'-matrix days, it was noticed Q] that a chiral Lagrangian used at tree level 
was a fleld theoretical realization of current algebra. Second, the spectacular experimental 
successes of the Standard Model together with its only partially unifled structure suggested 
that it be viewed as the low-energy remainder of a more fundamental theory [Q. Third, a 
new view of renormalization grew out of renormalization group ideas 0. By the end of the 
70 's it was clearly realized that the concept of effective Lagrangians provides the rationale 
behind the use of field theory in particle physics phenomenology P,01. It was in the 80's that 



^Currently, the only way to deal with processes with momenta Q ~ Mqcd is by phenomenological 
models. Those (few) models which include QCD symmetries correctly can be roughly thought of 
as EFTs in which an infinite number of operators is assumed to be correlated and their parameters 
functions of the few model parameters. Both model- independence and a systematic expansion are 
lost. Truly extending the EFT approach to such "high" momenta presents formidable obstacles. 
First, the number of relevant degrees of freedom increases tremendously. Second, one would need 
to find a new small expansion parameter. I do not attempt to review efforts in this direction here, 
as they have had no clear successes in nuclear physics yet. 



a detailed, systematic analysis of low-energy hadronic physics based on EFT began [§,3, 
but even then, no non-perturbative effects were dealt with. 

Here I will consider mainly the description of nuclear systems per se for Q < Mqcd, a 
program that was initiated in Refs. |ll0|-|l8|. This early work has shown that many of the 
observed but otherwise mysterious features of nuclear physics can be naturally understood 
in this framework |jT5|Jl^. More recently, several aspects of this program have received 
vigorous attention — for a sample, see Ref. [^ — and a reasonable body of successes has 
been amassed. A review seems timely. However, as many review writers before, I am 
compelled to apologize to those who feel that their work is underrepresented here. Spacetime 
constraints — the editor's patience notwithstanding — have biased me towards the ideas that 
I understand better, usually (but not always) my own. For a different perspective on various 
subsets of the topics I cover here, the reader is invited to consult Ref. pT| . 

After an introduction to the general ideas of the EFT method in Sect. |I[ I briefly present 
some xPT highlights in Sect. |TT1[ Sect. |V| on the two-body system is perhaps the central 
one. I consider first the very-low-momentum EFT as its simplicity makes it worth studying; 
then I delve into pions at low and moderate momentum. Power counting is the theme that 
has been played but whose scope is not yet completely understood. This sets the stage for 
the application to few-body systems in Sect. |V], where remarkable universality emerges. 
The question we explore is what are the new ingredients brought in by new bodies. Several 
successful applications to processes involving external probes are mentioned in Sect. 0, 



and the first, exploratory studies of the many-body system are in Sect. [VI I| . The balance is 
in Sect. |VlTTi 



II. EFT 



Effective field theories exploit the existence of scales in a system. A very brief introduc- 
tion to the main EFT ideas is presented Subsect. |II A| . In Subsect. [II B| we identify the 
relevant degrees of freedom whose propagation should be accounted for explicitly in order 
to supply the correct low-energy analytic structure of physical amplitudes. The symmetries 
which constrain the low-energy S*- matrix elements are discussed in Subsect. [II C[ . The 5*- 
matrix can then be generated by the most general Lagrangian involving the relevant degrees 
of freedom and symmetries. In Subsect. \\.\.D\ one invokes the concept of "natural" size for 



effective dimensionful parameters in terms of the scale M of the underlying dynamics, in 
order to establish an ordering of all possible contributions to processes at momenta of order 
Q -C M in powers of the small parameter Q/M. Some features of the method are stressed 
in Subsect. [HE. 



A. What is effective? 



What is an effective theory? In order to have a rough idea — a fuller answer can be 
found in several reviews p2|-[26[| — consider a field theory (which I call the "underlying" 
theory) given by some Lagrangian Cund written in terms of some ("elementary") fields \&; 



and suppose that this theory adequately describes experiments carried out over a certain 
range of energy 0. This means that S'-matrix elements can be obtained from the path integral 

Z = Iv^e'^'^--^^^^. (4) 

We are interested in how such a theory looks like at scales smaller than some scale A within 
its range of validity. Excitations with momenta larger than A cannot be produced directly, 
so we can profit from splitting the fields "^ in two components "^h ("fast") and \&/ ("slow"), 
according to whether their momenta are greater or smaller than A. Integrating over \E'/i, 

where the effective Lagrangian >Ce// is given, in D spacetime dimensions, by 

J d^'xC^jfi^i) = -zln|pv^,e^/^"-(*-*') = Jd^xJ2g,{A)0,{^i). (6) 

The operators Oi may be very complicated but they involve only the slowly- varying fields 
"^i; the fields \E'/i are said to have been "integrated out". The Oj's have two important prop- 
erties. First, they are local in the sense that they involve only fields at the same spacetime 
point. Yet, they contain arbitrary number of derivatives of such fields. Indeed, accord- 
ing to the uncertainty principle particles with momentum Q ^ A can only probe distances 
l/Q ^ l/A and therefore can sense configurations with support smaller than 1/A only in an 
average sense. Second, as long as the splitting of fields is done carefully, the Oj's transform 
under various groups according to the underlying theory. If the underlying Lagrangian is 
symmetric under a transformation that is not anomalous, then so is the effective Lagrangian, 
although the symmetry might be realized non-linearly if it is spontaneously broken. If the 
symmetry is explicitly broken either at a classical or a quantum level, operators will appear 
at low energies that incorporate the breaking accordingly. Note, however, that the splitting 
of fields does require some care. The Lagrangian Ceff typically contains interactions involv- 
ing the time derivative of "^i, "^i, to powers greater than 1. In this case there is a non-trivial 
relation between "^i and its conjugate momentum. As a result the interaction Hamiltonian is 
not simply minus the interaction Lagrangian, but contains additional, Lorentz non-covariant 
terms. It is no problem to include these interactions, as other non-covariant pieces arise in 
covariant perturbation theory from contractions involving derivatives, and in time-ordered 
formalism from its inherent non-covariance. One can show explicitly |^ that the sum of 
diagrams contributing to any given process is indeed covariant. Alternatively, one can add 
terms AC directly to C^ff- For example, for real scalar \E'/, if 

C^ff = hjAi^i)^i + . . . , (7) 



then [EB 



•f Most likely the underlying theory is itself an effective theory to an yet more fundamental theory. 



AC='-5''iO)Tr[ln{A-\^i))] + ... (8) 

This Lagrangian is manifestly covariant, but the new non-hnear piece generates self- 
interactions. These are certainly ill-defined as they stand due to the S^{0), but other terms 
in the unrenormalized Lagrangian are equally ill-defined. The EFT is made well-defined 
through regularization. S^{0) vanishes in dimensional regularization but not when a mass 
cutoff is used; the In term is necessary |^| in the renormalization program to get rid of loop 



contributions that could otherwise spoil the assumed symmetries. 

As for the coefficients Qi, they carry information about details of the dynamics: they are 
functions of the parameters of the underlying theory, but also depend explicitly on the cutoff 
A — which lends them the status of "running coupling constants" . In diagrams, a change in 
A amounts to a reshuffle between contributions from vertices and from the ultraviolet region 
of loop integration. An obvious requirement is that low-energy observables be insensitive to 
changes in A, which places constraints on the A dependence of the QiS. The equations that 
govern this dependence are the renormalization group equations. 

The notion of an effective Lagrangian is particularly useful when the underlying theory 
has (at least) one characteristic mass scale M. In this case, the effective degrees of freedom 
for A < M expressed by the fields "^i are in general significantly different from the origi- 
nal degrees represented by the fields "$. It frequently happens that that these low-energy 
phenomena are "collective excitations" as viewed in terms of \I' [^ . The reformulation of 
the theory in terms of "^i ends up selecting those effects of the underlying theory that are 
most important at low energies. The simplest case is the one where M is just the mass of 
a physical particle. Production and decay of this particle involve large momenta and do 
not concern the EFT. Effects of virtual exchange (of the particle, or of particle-antiparticle 
pairs) are of short range and thus only included indirectly in the Qi's [^]. If, furthermore, 
the particle is not stable in the context of the underlying theory, then it does not (explicitly) 
appear at all at low energies; we need associate no field ^^ to it. Another common case is 
that where M is the scale associated with some (elementary or composed) scalar field ac- 
quiring a non-vanishing vacuum expectation value and breaking a (possibly approximate) 
continuous global internal symmetry group, resulting in the appearance in the spectrum of 
either a (possibly nearly) massless spin- zero particle |^ — a (pseudo)Goldstone boson — or 



a massive, longitudinal component of a vector boson |^. A third case is more complicated: 
M is a scale where some wild non-perturbative phenomena dominate and force asymptotic 
states into singlet representations of a local group, an example being confinement in a non- 
abelian gauge theory with the appropriate matter content. It seems adequate in this case to 
restrict the effective Lagrangian to fields that are singlets under the gauge group. There are 



other possibilities (like the existence of a Fermi surface p3| , |2^ ), or combinations of various 



cases, but the important point is, the mass M provides a measure of the strength of the 
high-energy effects that appear in the coefficients ^^j's. This eventually paves the way to a 
systematic expansion in powers of Q/M. 

If we know and can solve the underlying theory, then we can calculate the Qi^s and obtain 
the complete form of Ceff to the accuracy of the approximation involved in the solution. But 
it frequently happens that we cannot solve the underlying theory (as it is the case of nuclear 
physics), or that we do not know such theory at all (as it is for the electroweak theory). Are 
we then justified in using an effective theory? A positive answer is given by a "theorem" 



due to Weinberg 0: in conciliating quantum mechanics and Lorentz invariance in a way 
consistent with unitarity, analyticity and cluster decomposition, we are led to quantum field 
theory; the most general Lagrangian with some assumed symmetries will then produce the 
most general S'-matrix incorporating those general principles and symmetries, without any 
other physical content. This is a statement based on our own currently accessible energy 
experience: it has only been proved in the particular case of a scalar field with Z2 symmetry 
in Euclidean space |]5^, but no counter-examples are known. 

EFT is therefore modern S'-matrix theory. All that any theory can do is to relate 
observables. Field theory based on a Lagrangian is the simplest way to do this given some 
degrees of freedom and symmetries. The theory is of course regulated by some sort of 
procedure implementing the cutoff A < M that isolates the neglected degrees of freedom, 
so diagrams involve only momenta Q < K. Renormalization is carried out as usual, but its 
aim is not to eliminate infinities: because the EFT is not valid at high energies, we never 
have to take the A — >■ 00 limit anyway |]23|. Rather, the objective is to remove details of 
the cutoff procedure from physical amplitudes, and to relate Lagrangian parameters and 
observable quantities |. A redefinition of fields in general changes all amplitudes off shell, 
where they cannot be measured. Only observables, which appear in on-shell amplitudes, are 
unambiguous in the EFT. 

B. Not too much freedom 

The first element of an EFT consists of the identification of relevant low-energy degrees 
of freedom. Here we will be considering only EFTs for momenta Q <^ Mqqd ~ 1 GeV. 
This severely restricts the number of degrees of freedom that need to be included explicitly: 
at such low momenta we cannot probe quark degrees directly, and the number of hadronic 
states needed is small. 

The nucleon mass rriN- is not light compared to Mqcd, but to a very good approximation 
nucleons cannot decay into lighter states of high momenta. Although we can integrate 
out nucleon-antinucleon pairs that pop out of the vacuum, a number of nucleons exist in 
asymptotic in- and out-states of nuclear systems. Nucleons of momenta Q <^ Mqcd are 
thus associated in the EFT with an explicit field N, and nucleon propagation is represented 
by a line that goes through diagrams. The EFT splits into sectors of definite nucleon 
number A: operators containing more than 2A nucleon fields do not contribute to processes 
involving only A incoming (and outgoing) nucleons. Because Q -C wijv, an incoming slow 
nucleon receives only little kicks from other particles. Nucleons are thus non-relativistic, 
characterized in first approximation by a non-relativistic dispersion law for the energy, E = 
p'^/2mN as function of the three- momentum p. Relativistic corrections are accounted for 
systematically in an expansion in Q/itim- 



§ This can only be done because we include all allowed interactions. This is in sharp contrast 
to what is usually done in models like the Nambu-Jona-Lasinio model of low-energy QCD; such 
models are frequently called "effective" , but they are not theories in the above sense. 



Which other hadronic states need to be included depends on what range of momenta we 
are interested in. Certainly all mesonic states with masses rrim ^ Mqcd can be accounted 
for indirectly, through their contribution to coefficients of local interactions. In complete 
analogy with the multipole expansion in classical electrodynamics, we can approximate the 
interactions among nucleons originating from the exchange of these mesons in a series of 
contact interactions with an increasing number of derivatives. In first approximation such 
contact interactions are momentum-independent, and corrections generate an expansion in 
Q/rrim- Likewise, baryonic states with masses uin* ^ ttin + Mqcd cannot be reached by 
hitting nucleons with small momenta. They can also be integrated out, generating local 
interactions in an expansion in Q/{m]^* — 771^). 

The most relevant state besides the nucleon is the lightest hadron, the pion of mass 
ttItt ~ Mnuc'i next is the delta isobar, as its mass splitting to the nucleon is not large, 
5m = m^ — mjsf ~ Mnuc- For Q ^ M„„c the pion and delta can be integrated out — unless 
we are considering external pions, in which case they could be treated as heavy particles 
as well. For Q ^ 711^^ the relativistic pion has to be included as an explicit field tt. For 
Q ^ 771/^ — 77iN a non-relativistic delta field A has to be considered as well. Because the 
two thresholds do not exactly coincide, it is possible — and possibly efficient — to integrate 
out the delta in a limited region of momenta above the pion threshold, while maintaining 
explicit pion fields. Coefficients of the corresponding effective Lagrangian that receive direct 
contributions from the delta will be relatively large, being suppressed only by powers of Q/5. 

C. Symmetries 

Symmetries play a fundamental role in EFTs because they restrict the possible form 
of interactions. The EFT corresponding to the standard electroweak theory at an energy 
scale of a few GeV involves only the lightest leptons and quarks, gluons, and the photon; 
weak gauge bosons, and heavy leptons and quarks can be integrated out in favor of non- 
renormalizable interactions. Here I am interested mostly in interactions involving the lightest 
u and d quarks. They can be arranged in a ffavor doublet q = (^). Denoting by G^ (A^) the 
gluon (photon) ffeld of strength G^^ (F'^v), the relevant pieces of the effective Lagrangian at 
this scale are 

C = -g(^ - igs(p- ieQ4)q - -("^m + md)qq + -(m^ - m„)gT3g 






-Tr[G^,Gn + ^^e^,p„Tr[G^-G'"'] - -F^.F^^^ + ... (9) 



Here rs is the usual Pauli matrix, Q = 1/6 + T3/2 is the quark charge matrix, and ". . ." 
denote non-renormalizable terms. Non-renormalizable interactions are suppressed by powers 
of the masses of the heavy particles that have been integrated out. I will neglect them in 
a (very good) first approximation. The leading Lagrangian (QCD + QED of quarks) has 
5 parameters: the gauge couplings for strong (qs) and electromagnetic (e) interactions, the 
masses of the up (m^) and the down (m^) quarks, and the strong CP parameter (6). The 
theta term is found to be unnaturally small (the so-called strong CP problem), and I will 
neglect it as well. 

9 



The set of remaining interactions in Eq. (^ has a number of approximate symmetries. 
They are clearly invariant under proper Lorentz transformations, parity, and time- reversal. 
Local changes of the quark-field phase that are proportional to unity or generators of the 
SU{3) color group leave the Lagrangian invariant. A global U{1) symmetry corresponds to 
baryon number, while another, axial U{1) is anomalous. In the limit where m„, rrid, and e 
are zero there is further invariance under a chiral symmetry SU{2) x SU{2), consisting of 
independent rotations with parameters Sl and Sr of the left (g^ = (1 + 75/2)^) and right 
{qr = (1 - 75/2)g) quark fields: 

5qL = -iSL ■ tqi, 5qR = -iSr ■ tqR, (10) 

where t = t/2 are the generators of SU{2) in terms of the Pauli matrices r, 

1 i 

iah = -T^ab + -^abctc- (H) 

Acting on quark bilinears, SU{2) x SU{2) ~ 5*0(4). The absence of degenerate parity 
doublets but presence of (approximate) isospin multiplets in the hadron spectrum indicates 
that chiral symmetry is broken by the vacuum down to the diagonal subgroup SU{2)v ~ 
5*0(3) of isospin. Mass terms and electromagnetic interactions are not invariant under chiral 
transformations: they transform as non-trivial 50(4) tensors and therefore break chiral 
symmetry explicitly. The operator qq breaks chiral symmetry as the fourth component of 
an 5*0(4) vector 

S={2qtj,tq,qq), (12) 

while —qr^q is the third component of another 50(4) vector, 

P=i-2qtq,qi^,q). (13) 

Photon exchange among quarks produces four-quark interactions off which we can read the 
way electromagnetic interactions break isospin. Defining the 50(4) antisymmetric tensor 

\ -qi-f%q J ' ^ ^ 

one finds that chiral symmetry is explicitly broken by terms proportional to T^^ and T^4T^^. 

The EFT has to produce the same low-energy 5'-matrix elements as the underlying theory. 
Although any choice of fields is equally valid, it is convenient to choose fields that transform 
under the above symmetries in the most transparent way. Although through field redefini- 
tions other forms of interactions can be obtained where symmetries arise from a conspiracy 
of several terms, the particular choices of fields described next are the most convenient, as 
the pattern of symmetries is respected term-by-term in the resulting Lagrangian. 

Spacetime symmetries. Parity and time-reversal can be implemented at low energies 
in the usual way, but proper Lorentz invariance deserves some care. Not surprisingly, the pion 
is represented by a relativistic pseudoscalar field with standard transformation properties, 
but heavy stable fields introduce a small complication. If I denote nucleons and deltas by 
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a generic field ip of mass m ~ M , tlien interactions will contain c^o'?/' ~ —imip + . . . and 
produce factors of m/M that need to be resummed before we can develop an expansion in 
Q/M. On the other hand, the large energy m is not really available for the dynamics, so we 
expect the result of the resummation to be well behaved. In fact, the result can be no other 
than a rearrangement of the non-relativistic expansion of the relativistic theory. 

The resummation can be done from the start, avoiding the relativistic theory altogether, 
in what is referred to as the heavy-particle formalism p4| , ^ . We extract the rest energy by 
a field redefinition, if) = e~^"^'"'^ilj in a frame where the heavy particle moves with velocity 
V. The most general Lagrangian involving ip will be free of explicit factors of m/M but 
will contain explicit factors of f^. Proper Lorentz invariance requires that the Lagrangian 
be invariant under a boost to a frame of different velocity w^^ = f^ + Q^/iTi, with q <^ m 
satisfying v -q = —q^ /2m. Fields in the two frames are related by a phase, besides the usual 
matrix D{q) that accounts for the appropriate representation of -^ under the Lorentz group: 
'4'w = e^'^'^D{q)ip^. Derivatives of ip transform in a more complicate way; as usual in such 
circumstances, it is profitable to introduce a covariant derivative 

such that d^^^ transform in the same way as ipy itself. Out of d^^ip and (d^^)^ = ip'^dj^ we 
can construct a covariant velocity and a symmetric derivative 

V. = ^(d,-dt), ^,^^(d, + dt), (16) 

which are useful building blocks, as in bilinears they produce Lorentz invariants. Spin can 
be dealt with in similar manner. If ip has spin 1/2, then at low energy only two degrees 
of freedom are important, which can be enforced by imposing that the spinor ip^ obeys 
(1— /j)il'v = 0. The covariant spin is 

S^ = '-l,a,,V^, (17) 

as indeed E ■ V = 0. If ^ has spin 3/2, then at low energy only four degrees of freedom are 
important; this can be enforced by demanding the field {ipv)fi (which behaves as the tensor 
product of a Dirac spinor and a four-vector) obey v ■ ^py = and '~i^o'^"v^{ipv)v = 0. To 
account for the vector character, the covariant spin-3/2 field has to be 

(^,)^ = i^,)^ - -v^d ■ ^,. (18) 

m, 

Lorentz invariance is ensured by proper contraction of the indices of these covariant 
objects ||3^. These invariants reduce in the rest frame of the heavy particle to a set of 5*0(3)- 



invariant operators, all with relative strengths determined by 1/m. The heavy-particle 
formalism is but a way to impose invariance under rotations plus slow velocity boosts; 
truncated at leading order in 1/m it is just familiar Galilean invariance. As an example 
we can look at kinetic terms for ip, i.e. terms containing only two heavy fields. Any such 
bilinear involving (}^ can be eliminated in favor of total derivatives, so the most general 
kinetic term is of the form ^^^(V^ — 1)"^,;, n > 0. By field redefinitions involving (V^ — 1) 
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we can eliminate all but ipliiv ■ d — d'^ /2m)ipv] by field redefinitions involving iv ■ d and d^ 
we can systematically eliminate powers oi v ■ d bigger than one. In the rest frame v = (1,0) 
we then arrive at the standard 

r = ,.tU + _Lv^ + _i_v* + ...),, (19) 



When two heavy fields of mass m and m + 6m, 6m -C M, are considered in the same EFT, 
it is simpler to redefine both with the same phase e*™^'^. In this case the kinetic terms of 
the heavier field are as in Eq. (^) with additional terms —ip'^6m,{l + . . .)?/'. 

A consequence is that we can alternatively work in the rest frame with heavy fields that 
provide definite representations of 5*0(3) ||15|. By field redefinitions we can systematically 
remove time- derivatives of the heavy fields from interactions. We take A^ to be a Pauli spinor 
that provides the 1/2 representation |(j of the spin generators, 

ajCTj = 6ij + ieijkCTk- (20) 

In the rest frame, ip = ([)N/\^. Likewise, the A provides the 3/2 representation of the spin 
generators. In the rest frame, i/j^ = (0,S'A), 

SiSf = -{26ij - ieijkCTk)- (21) 

o 

The 2x4 transition matrix S is used to couple A to A^ in spin-1 bilinears. 

Internal symmetries. Goldstone's theorem |31[] tells us that massless Goldstone bosons 
arise due to the assumed spontaneous breaking SO (4) — > 5*0(3), and are naturally identified 
as pions. Pions are then associated with broken generators of 50(4), and their fields live in 
the three-sphere 5'0(4)/50(3) ~ 5*^. We call the diameter of this sphere F^; it is directly 
related to the chiral symmetry breaking scale, ^xSB, but the precise factor can only be 
obtained from the (so far elusive) mechanism of dynamical symmetry breaking in QCD. 
It turns out that this diameter can be determined from pion decay, and is the pion decay 
constant, F^^ ~ 190MeV. If we embed S^ in the Euclidean E^ space, 5*0(4) transformations 
can be viewed as rotations in various E^ hyperplanes. For example, SU{2)v of isospin 
consists of rotations in planes orthogonal to the fourth axis, while axial SU{2)a are rotations 
through planes that contain the fourth axis. 

The sphere can be parametrized any way you want, say with four cartesian coordinates 
{(p, (y94 = a} subject to the constraint a"^ + cp"^ = -F^/4. The coordinates {(p, a} transform 
linearly under 5*0(4) rotations; for example, under an infinitesimal SU{2)y transformation 
with parameter Sy, 

6v(fi = By y< (p. (22) 

It is more convenient, however, to work with three unconstrained coordinates tt that can be 
directly associated with the pions. Any point on the sphere can be obtained by applying a 
four-rotation R{7v), i?(7r)i?^(7r) = 1, to the north pole (0, ^F.„): 

<^a(7r)=i?„4(7r)y. (23) 
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Different parainetrizations of R define different pion fields tt. The tt's just rotate, 

(5y7r = £v y< 77, (24) 

under SU{2)v, but transform non-linearly under SU{2)a_. From a particular choice of -R(7r) 
one can read the non-linear transformation of tt. The transformation law for (9^7r is invari- 
ably complicated. 

The general theory of non-linear realizations of a group is discussed in Ref. I^Gj. The 
main idea is to construct covariant objects that transform under the broken generators as 
under the unbroken subgroup, but with a field-dependent — thus local — parameter. The 
transformation we want to imitate is the rotation ( ^2]) and the appropriate field-dependent 
parameter for an axial transformation with parameter e^ is Ea x tt/F^. We demand that 
the covariant derivative of the pion field 

D, = ^{l + 0{77yF:)) (25) 

transform as an isospin-1 object. 



TT 



SvDi, = sv X Df,, SaD,, = [sa x ^r) x D^. (26) 

Given the local look of 5aD^, it is not surprising that the covariant derivative of D^ is 
formed from a vector field 

^. = |^^x9,7r(l + 0(7rVF,2)) (27) 

as 

V,,D, = d^D, + iE^ X D,. (28) 

The pion fields being associated with broken generators, one can amass them in a matrix 
M(7r ■ t) = a(7r^) + ^6(77^)77 ■ t — with a(7r^) and b{n'^) two functions that depend on -R(7r) — 
defined so that 

i 1 

D^-t = --(u^d^u-ud^u^), E^-t = -{u^d^u + ud^u^). (29) 



For example, one can use stereographic coordinates |TD| tt = 2c^/(l + 2a /F^), which 
correspond to 

R[-]=\^" on-l^/^^ n-1^ \^\ r (30) 

where 

-,r2 

(31) 

In this case, the pion covariant derivative and the vector field take particularly simple forms 
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-2D-i"^^- 


2D-^^ 


2D-i^ D 


-1 /1 7r2 


o-^^.t- 





D, = D-^^-f^ E, = ^.xD,. (32) 

-^ TV -^ 7T 

Other common choices of pion field tt are the "sigma-model" parametrization 



u' 



47r2 Ai .. 

+ — TT ■ t (33) 



\ F^ F. 



and the "exponential" parametrization 

u = e^^ -^ (34) 

Fields that couple to pions can be introduced easily once their isospin character is known. 
It simplifies to work with fields that transform under broken rotations with the same field- 
dependent parameter as before. In the case of a (2n + 2)-component field %jj^"'^ providing a 
representation f ("■+^/2) of isospin n + 1/2, 

Its covariant derivative is then 

^M^^"^ = (5m + *^"^'^'^ ■ ^m)^^"^ • (36) 

A nucleon field A^ is a two-component isospinor and Eqs. (|35[) and (|36|) hold for i})^^^ = N 
and t^^'"^^ = t. A delta field is a four- component spinor A and the same formulas apply for 
tp^^^ = A and t^^'"^' . Isospin-1 bilinears of A^ and A can be formed using a 2 x 4 transition 
matrix T that satisfies 

TaTf^ = -{Sab - ie abate)- (37) 



The introduction of fields that transform non-linearly makes it easy to construct chiral 
invariant objects: any isoscalar built out of the fields D^ and %p, and their covariant deriva- 
tives will automatically be invariant under the whole SU{2) x SU{2) group. Operators 
Ta/3...['^', D^, i'] that break chiral symmetry as tensor products of the tensors (0), (p^, and 



(|l^) must also be constructed so that we can reproduce the correct symmetry pattern of 
low-energy S'-matrix elements. Analogously to Eq. (P5D, where an SO {4) vector involving 
the pion field tt is constructed by applying a chiral rotation on another that does not, we 
just need to search for those tensors Ta^...[0; Z?^, ip] written in terms covariant objects only, 
then rotate them with R{tt): 

Tap...[7z;D^,i;]= ^ i?«„/[7r]i?^/3/[7r] . . . Ta'p'...[0; D^,ip]. (38) 

a'/3'... 

For example, the simplest vector we can construct out of numbers only is 5'(i)[0; 0, 0] = (0, 1) 
so from (^) and (|30D , the fourth component of 

SiiM;O,0] = (^2D-'^^,1-2D-'^^ (39) 
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breaks chiral symmetry in exactly the same way as qq. By considering other vectors and 
tensors we can construct the infinite number of operators that break chiral symmetry as 
quark masses and electromagnetic interactions. 

Because the EFT is formulated in terms of the lowest-lying color singlets, color gauge 
invariance plays no direct role. Soft photons (those with momenta < M) are not inte- 
grated out of the EFT and couple to pions, nucleons, and deltas in the most general way 
that respects U{l)em gauge invariance. The photon field A^ enters through (i) "minimal 
substitution" in covariant derivatives, 

^^,^a -> d^,^a - eA^egabVTfe, d^i)'^'''^ -^ d^ip - ieAf,Q^ip^''\ (40) 

where Q^p is the charge matrix of the field ip^'^^ — for example, Q^ = (1 + ''"3)/2; and (ii) 
gauge invariant objects constructed out of the field strength F^^. There is, however, one 
subtlety regarding the role of gauge invariance: chiral symmetry is broken explicitly by the 
electromagnetic interaction because there is an anomaly in the third (isospin) component 
of the axial current. Terms have to be included in the effective Lagrangian in such a way 
that under a chiral transformation they reproduce the anomaly. This abnormal intrinsic 
parity sector is reviewed in Ref . |^ . Finally, further breaking of the above symmetries by 
the theta term pSf and non-renormalizable interactions |39| in the Lagrangian (P) can be 



included as above, mutatis mutandis. 

D. The power of counting 

Although symmetries restrict the number of interactions in the EFT, we still have to deal 
with an infinite number of them. The only hope of any predictive power relies on finding 
some ordering of contributions. Since it is not likely that the infinite number of coupling 
constants will be all small enough for a coupling-constant expansion to converge, we turn to 
the only obvious small quantity, the typical momentum Q itself. 

Powers of Q of any particular Feynman diagram can be counted in ways analogous to 
the one used in finding the superficial degree of divergence (see, e.g., Ref. [0). Each space 
derivative in an interaction produces a three-momentum in a vertex and therefore counts as 
Q. Likewise, each time derivative of a light field brings a Q as well. In any loop, integration 
over the zeroth component of the four-momentum will generate two types of poles according 
to the scales appearing in the propagators: (i) standard poles at ~ Q corresponding to 
external three-momenta and to the masses of the lightest excitations; (ii) shallow poles at 
~ Q"^ /2m]^ corresponding to initial energies. 

Processes that involve at most one heavy particle line are the simplest because we can 
always close the contour avoiding shallow poles. As a consequence, each four-momentum 
integration brings a factor Q^. A nucleon propagator is Q~^ from the first term in the kinetic 
Lagrangian (0), the other terms being treated as corrections: the nucleon is a nearly static 
source propagating forward in time, as can be seen from the Fourier transform 9{t)5^{r) + . . . 
of its propagator. If Q is sufficiently high so that pions and deltas are included explicitly, a 
pion propagator is Q~^, and a delta propagator is Q~^. 

Processes that involve more than one stable heavy particle present a complication: a 
failure of perturbation theory that can potentially lead to bound states [^]. To see this, 
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(a) 



(b) 



FIG. 1. (a) A two-nucleon intermediate state between two unspecified interactions, associ- 
ated with a loop four-momentum /; (b) the fuh two-nucleon amplitude T^n as an iteration of a 
two-nucleon potential Vnn- 

consider the scattering of two heavy particles with center-of-mass momentum p and total 
energy E = k'^/m + . . .. In-between two interactions, the heavy particles propagate forward 
in time. Such a piece of an arbitrary diagram is shown in Fig. |l|(a). Neglecting smaller 
relativistic corrections, it will contribute 



rf/o 



1 



1 



2n p^E 



2m K 



+ %t -P + f 



2-mjv 



+ ie 



(41) 



to the result of the diagram. The problem is that this integral is infrared enhanced, for 
it picks a large residue from one of the shallow poles. A naive momentum counting would 
have suggested a contribution of 0(1/(5), yet we end up with something of 0{mi^ /Q"^). 
This is quite a general phenomenon: such an enhancement of 0{M/Q) appears whenever 
two or more particles of mass 0{M) but momenta 0{Q) collide. As a consequence, for 
those contributions that come from shallow poles each four-momentum integration effectively 
brings a factor Q^. A nucleon propagator is counted as Q~^, and the first two terms in the 
kinetic Lagrangian (|1^) are important. A pion propagator still counts as Q~^, but the pion 
can be taken in first approximation as static, and it is sometimes referred to as a "potential" 
pion. Contributions that come from standard poles naively scale as in processes with none 
or one heavy particle. Pions there are non-static or "radiative". 

After integration over energies, each graph in covariant perturbation theory produces 
a number of terms. One can associate these with graphs in time-ordered perturbation 
theory, where loops represent three-momentum integrations and internal lines stand for 
denominators which contain energy differences between initial and intermediate states. In 
general, we can distinguish two classes of time-ordered diagrams. (Ir)reducible diagrams are 
those that can(not) be split in two pieces if in an intermediate state we cut all — and only — 
the lines corresponding to initial or final particles; their denominators will (not) contain 
energy differences of 0(Q^/mjv). 

In a system of heavy particles, we call the sum of irreducible diagrams the potential 
V. Note that the potential being a set of subgraphs, it can be defined in alternative ways. 
All definitions that exclude the infrared enhancement contributions but differ by a smaller 
amount are equally good, as long as no double-counting or omissions are made. The impor- 
tant point is that the only scale appearing explicitly in the potential is Q, so that the power 
counting proceeds as in the case of diagrams with at most one heavy particle. 
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Reducible diagrams can be obtained by sewing together irreducible diagrams. The full 
amplitude T is a sum of the potential and its iterations; schematically, with Gq denoting 
the intermediate state containing only initial or final particles, 

T = -{V + VGoV + VGoVGoV + ...) = -V + VGqT. (42) 

This is just the Lippmann-Schwinger equation, which is formally equivalent to the 
Schrodinger equation with the potential V. Fig. |I](b) illustrates the case of two nucleons, 
which formally evolve according to the familiar non-relativistic Schrodinger propagator 

Go(/;A:) = -^^-^^, Go(r; A:) = -m^. (43) 

Relativistic corrections are accounted for by a two-legged vertex 

/4 - A-4 
u{l;k) = -—^ + ..., (44) 

and can be lumped in the potential. 

With these ingredients one can write the contribution of any diagram to the amplitude 

as 

TocQ'J^iQ/A), (45) 

where A is a renormalization scale, and i/ is a counting index. The exact formula for u 
in terms of the elements of a diagram depends on the Q we are interested in. In addi- 
tion, electromagnetic interactions can be considered through a simultaneous expansion in 
a = e^/47r. We will discuss the power counting in more detail below. 

Counting powers of Q is not in itself sufficient for an ordering of interactions. We need 
to find which other scales accompany Q. As we have mentioned in the introduction, there 
are at least three physical scales in the problem: K, Mnuc, and Mqcd- The largest possible 
ratio is Q/K. From dimensional analysis alone, then, we expect that for Q ^ K amplitudes 
for scattering of free nucleons will be perturbative in QfR. Such systems are labeled dilute. 
As Q increases, we have to perform resummations: first summing all orders in Q/K but 
expanding in Q/Mnuc] then summing all orders in Q /Mnuc but expanding in Q/Mqcd] and 
so on. As we resum certain class of diagrams, bound-states may appear. 

Given that we do not yet possess a full solution of QCD, it is obvious that some assump- 
tions have to be made about the parameters of the EFT. The minimal assumption is that 
of naturalness PTl . Naturalness means simply that for each operator Oi of dimension 5i in 



Eq. (H), the parameter gi{k) will be of form Ci{k/M)/M^' with Ci(l) = 0(1) Q, unless there 
exists a symmetry that forces Cj to be small. As we are going to see, the presence in the 
nuclear problem of scales lighter than Mqcd requires slightly stronger assumptions about 
the scaling of the EFT parameters. These new assumptions are still fairly general and are 



**Using a mass cutoff A <C Af introduces a large dimensionless ratio M/A and can unnecessarily 
complicate the power counting. 
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motivated by qualitative features of nuclear phenomenology; they are much less restrictive 
than usual model assumptions. In any case, these assumptions are necessary for the a pri- 
ori ordering of contributions, but they are subject to test a posteriori. If a parameter is 
found much larger than assumed initially, the EFT can be improved to better account for 
its effects. 



E. Propaganda 

The EFT approach is: 

(z) systematic. Once a power counting is established, there is an expansion in a small 
parameter oc Q. 

(a) consistent with QCD. The only (very important) QCD inputs are confinement, (chi- 
ral, ...) symmetries, and some version of naturalness. QCD can be represented by a point in 
the space of renormalized parameters (at some renormalization scale) of the effective theory. 
An explicit solution of QCD (such as from simulations on a lattice) would provide knowledge 
of the exact position of this point, and then the EFT would be completely predictive. Until 
such a solution is found — or as a test of QCD after it is found — we can recourse to fitting 
of low-energy experiments in order to determine the region in parameter space allowed on 
phenomenological grounds. Even in this case the theory is predictive, because to any given 
order the space of parameters is finite-dimensional. After a finite number of experimental 
results are used, an infinite number of others can be predicted up to an accuracy depending 
on the order of the expansion. In practice, because the number of parameters grows rapidly 
with the order, model-dependent estimates of parameters based on specific dynamical ideas 



— such as saturation by tree-level resonance exchange [|^ — are sometimes used. 

{Hi) a consistent quantum field theory. Renormalization can be carried out to relate 
parameters appearing in different processes and to remove cutoff dependence at each order 
in the expansion. As a consequence, it is applicable in principle to all low-energy phenomena. 
The literature on applications to mesons and one-nucleon processes is vast; after a lightning 
review I will show how we can also gain insight into nuclear physics per se. 



III. ONE- OR NO-NUCLEON PROCESSES 

The EFT at momenta Q <^ Mnuc contains only nucleons, pions treated as heavy particles, 
and photons. All interactions are of contact type, and form an expansion in Q/Mnuc- The 
relevant symmetries are Lorentz and f/(l)em local invariance. Processes involving two or 
more pions, or one nucleon and at least one pion are of the type considered in the next 
sections. Processes involving one pion or one nucleon and a number of photons give rise to 
the various well-known electromagnetic parameters. Particle form factors take a multipole 
form, that is, an expansion in charge and anomalous magnetic moment, radii, etc. Compton 
scattering has a Thomson-seagull interaction and a magnetic moment piece, followed by 
polarizabilities, and so on. The analytic nature of these expansions in Q/Mnuc and the 
limited symmetry constraints make this EFT not very interesting. We expect, of course, 
that a sizable fraction of the value of the parameters comes from the lightest state that was 



integrated out, the pion. As Q is increased, the approximation of zero-range propagation 
for pions becomes less and less reliable. 

At momenta Q ~ m^, pions have to be included explicitly in the theory. Because 
numerically 5m is ~ 2?7i^, convergence of the EFT for Q ~ Mnuc should be optimized by the 
concomitant inclusion of an explicit delta degree of freedom. Now, we are in luck because 
pion interactions are not arbitrary. Once explicit pion fields are considered, approximate 
chiral symmetry imposes important restrictions on the way pions couple among themselves 
and to other degrees of freedom. The delta, too, can be included without too much hassle, 
since at these momenta it is, like the nucleon, a non-relativistic object, to which we can 
apply the heavy-particle formalism. In the preceding section all the ingredients to construct 
this EFT were presented. 

The power counting ( ^Sj ) presents no problems as the only explicit scales are Mnuc and 
Q, contributions from momenta of order Mqqd being lumped in the vertices. Assuming 
Q ~ Mnuci it follows that for a diagram involving A = 0, 1 nucleon and any number of pions 



2-A + 2L + Yy^^u ^i = d, + f,/2 - 2, 



(46) 



where L is the number of loops and Vi is the number of vertices of type z, which contains di 
derivatives or powers of m^r and fi fermion fields P,p!0| . This formula is important because 
chiral symmetry places a lower bound on the interaction index Aj > 0. Since L is bounded 
from below (L > 0), z/ > Vmin = 2 ~ A for strong interactions. An expansion in Q/Mqcd 
results. It starts at z/ = I'mtn with tree {L = 0) diagrams built out of vertices of index 
(X^V^Aj = 0), then proceeds at z/ = Umin + 1 with further tree diagrams, now with one 
vertex of index 1, the remaining having index {J^Vi^i = !)• These first two orders are 
equivalent to the current algebra of the 60's, but now unitarity corrections can be accounted 
for systematically. At z/ = z/^m + 2, for example, besides tree diagrams with one index- 
2 interaction or two index-1 interactions {J^^i^i = 2), there are also one-loop (L = 1) 
diagrams built out of index-0 vertices {J2 VJAj = 0). This is generalized to higher orders in 
obvious fashion. We discuss isospin breaking soon. 

It is convenient then to split the chiral Lagrangian in pieces C^^^ labeled by the index 
A. For example, the lower-order Lagrangians are |1T5|JT6| , |43| , ^ 



£(0) = i£)2 _ _^2^-i7r2 + nHVqN - 2gAt ■ a ■ DN 
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(48) 



Here qa, h-A = 0(1) and Bi = 0{1/Mqcd) are undetermined constants, to be obtained either 
by solving QCD or by fitting data; ". . . " stand for other terms involving the delta isobar 
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which do not appear exphcitly below. Higher-index interactions can be constructed similarly 
Note that pion self-interactions appear only at even index because of Lorentz invariance. 



The most general isospin-violating Lagrangian |l5| , p!8| is constructed out of operators 
that break isospin in three different ways. First, there are gauge invariant couplings to soft 
photons. For example, from minimal substitution in the index-0 Lagrangian, 



Ci-;'^ = -eA^^in X d,7z) 



-6^2. 2 

2 ^ 



VTn 



-eAoN^l + ts)N + '^A ■ N^a{t x 7z)sN + 

2 r^r 



(49) 



Second, there are operators that transform as products of the quark-mass-difference term; 
they are proportional to ewi^, where e = rrid — mu/{md + rriu) ~ 1/3 |^^. For example. 



4!^ = SmN{-N%N + 2D-'^ ■ NHN) + . . . , 



(50) 



C?l = ^D-'Smlnl + (3,{D, 



2D- 



F3 



-TTsTT ■ D) ■ N^aN + 



(51) 



0{6ml/M^ 



QCDJ 



Third, 



with 6mN = 0{eml/MQCD), Sml = 0{e'^ml/M^c^)^ and A 

there are operators that transform as four-quark interactions generated by hard-photon 

exchange, which are proportional to a. For example. 



C 



(-2) 
hp 



--D~'6ml(7z^ - nl] 



(52) 



'-'hp 



6m 



N 



2D 



-1 



NhsN + ——N\7Z% - TTgTr ■ t)N 



p2 



+ 



(53) 



with 6ml oc aMQ(^j^ and 6m.N = 0{6m.l/MQCD)- 

The observed smallness of isospin violation at low energies is not immediately evident 
within QCD, as chiral symmetry breaking effects are not negligible, and e is not particularly 
small. Power counting offers an explanation for this phenomenon [^,|l^. The first task is 
to compare the different towers of isospin-violating operators. This can be accomplished by 
noticing that photon loops (containing two electromagnetic interactions) typically produce 



factors of a/ir, which is numerically ~ e{mT^/mpY . This immediately suggests that C 
£ll+3) and C^rJ 



(n) 
hp 



■^qm 



sp 



(n+-) 

Cqm ^ . It can easily be verified that this estimate produces the right 
order of magnitude for the various contributions to the pion and nucleon mass differences, 
Am^ = 5m^ + (5m^ + . . . and Am at = 6m,N + 6m]^ + . . ., respectively. One now can see that, 
because there is no lowest-order isospin-violating interactions, isospin-violating quantities are 
typically not only 0{e) smaller than chiral- violating quantities, but are in general further 
suppressed by factors of Q /Mqcd'- isospin is an accidental symmetry [ll5| , p!8| . For simplicity. 



when considering only (a fixed number of) external photons, I will count a power of e as a 
contribution of 1 to the index A in Eq. (|46|). 
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Because 5m ~ Mnuc, delta effects appear for Q ~ Mnuc typically at the same order as 
nucleon contributions f^. At momenta Q < 6m we can integrate out the delta, producing 
an EFT with nucleons and pions only. This entails significant reduction in the number of 
interactions, as the delta terms disappear and the remaining interactions among nucleons 
and pions have the same structure as those in the C^^'^s above. Only their coefficients 
change, and now contain 0(l/5m) terms. For example, C^^^ has the same pion-nucleon 
seagulls with coefficients traditionally called Cj in the literature [H3| and related to the Bi 
via 



5, 1 / o2 N ^2 



^^ = T+---' '' = -ir'-2^r^s^^-- 



Bi hi 1 ( \ \ h\ , , 



4 9 (5m 4 V rrifq ) 18 5m 

This machinery can be applied to soft processes involving any number of pions and 
photons, and at most one nucleon. xPT is particularly relevant for those processes where 
the pion plays a dominant role. Here I briefly present some of the highlights that are 



directly relevant for what follows. For details on pion processes see, for example, Ref. [46 



one-nucleon processes in the deltaless theory have been extensively reviewed in Ref. [^ 



and explicit delta effects are discussed in Ref. [^ . 

Pion-nucleon scattering is perhaps the most important process where the chiral La- 
grangian can be applied, as it can be used to test the fermionic sector of the theory and de- 
termine some of its parameters. Contributions start at z/ = z/^m = 1 from Born graphs with 
intermediate nucleons and deltas, and from the ttttNN interaction ( "Weinberg- Tomozawa 
seagull" ) generated by the chiral covariant derivative in the kinetic nucleon term of Eq. (^7|). 
At z/ = Vrnin + 1 the Sj-seaguU diagrams contribute at tree level. Dropping nucleonic Born 



terms, the vr^A^ — ^ vr^A^ amplitude is traditionally written [^] as 

a + h q ■ q' + c (g ^ + q'"^) + e ojuj' — d eabc^c cr ■ qx q' + . . . (55) 



tab oc 6, 



ab 



where q (g ') is the incoming (outgoing) pion momentum and uo {uo') is the initial (final) pion 
energy. The a, ..., e coefficients are given in terms of the q's as 

IGmld 8c3 4c4 8(c2 + C3) 

TV TT TT TT 

At z/ = Umin + 2 a host of terms enter, including loop diagrams and new parameters. Very 
close to threshold, where only S waves matter, chiral symmetry predicts an isovector scat- 
tering length b^^^ ~ m^T^/F^ from the Weinberg- Tomozawa seagull [^ and a smaller isoscalar 



^TSome people reserve the name xPT to the EFT of the world where Q/5m ~ m-,^jbm <C 1, and 
refer to the EFT of a world where Q/5m ~ m-^jbm ~ 1 as the "small scale expansion". I find it 
simpler to refer to them as deltaless and deltaful xPT, respectively. 
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scattering length b^^^ ~ ml/F^MqcD, related to the sigma-term, from the subleading Cj seag- 
ulls. The most extensive calculations have been performed to z/ = Urnin + 2 in deltaless xPT. 
Different pieces of vrA^ scattering data have been used to determine the chiral coefficients. 
In lowest order, qa — 1.25; in sub-subleading order chiral symmetry breaking corrections 
appear, and if they are absorbed in qa, then qa — 1.32 satisfies the Goldberger-Treiman 
relation without discrepancy. In Table | we list some of the determinations of the param- 
eters Cj. Earlier fits ||43J5O|-02i were made to different sets of threshold and sub-threshold 



parameters obtained from dispersion analyses of older data. Newer fits [^ were made to 
different phase-shift analyses (PSAs), the last two in Table | including the more modern 
meson-factory data. The 0{Q^) determinations are consistent with each other when their 
error bars (not shown) are considered, except for ci, which refiects the higher value for the 
sigma-term in the newer PSAs. Note that the coefficients C2, C3, and C4, which receive con- 
tributions from the A at tree level, are larger than ci, as expected. There are not yet very 
reliable determinations of Ha from vrA^ scattering in the deltaful theory. The large- A/'c-limit 
value, Ha = Sqa/V^ — 2.7, gives a reasonable delta decay width at tree level. 

The only cases where isospin violation can be large are instances where the leading 
isospin-conserving interactions do not contribute. In such cases we can encounter an isospin 
violation of size comparable to e. One example is S'-wave ttN scattering [^,^,0,Q. The 
leading isospin-violating term involving a nucleon is, according to the above discussion, 
given by Eq. ([50|). It determines not only the most important quark- mass contribution 
SniN to the nucleon mass splitting, but also isospin-breaking interactions of an even number 
of pions with a nucleon. It contributes ~ 6mN/F^ to the pion-nucleon scattering lengths. 
This gives a small contribution ~ 6171^/171.,, ~ em-^/MQCD ~ 5% to processes that involve 
at most one 7r°, as they receive a leading-order contribution from the Weinberg- Tomozawa 
seagull. But for vr'^A^ — ^ tt^N, which starts at subleading order with the isoscalar seagulls, 
isospin violation can be as large as e ~ 30%, thus revealing the full isospin breaking in the 
QCD Lagrangian! Unfortunately, this one case of large isospin violation cannot be measured 
easily. 

Pion photoproduction at threshold, where the pion is mostly in an S wave, has received 
a lot of attention lately from both theoretical and experimental standpoints. This process 
has been studied near threshold up to z/ = Umin + 3 in the deltaless theory [HH| . At leading 



TABLE I. Subleading xPT coefficients in GeV ^ from several recent fits. 
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1.53 
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5.25 


4.12 


5.29 


3.63 


5.54 


3.25 


5.93 


3.44 


6.00 


3.52 


6.19 


3.51 
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order, u = Vmin = 1, there exist only contributions to charged channels from the fourth term 
in Eq. (^91) ( "Kroll-Ruderman term") and from the pion pole; we thus expect the charged 
channels to be larger. At z/ = i^rnin + 1, minimal substitution on the recoil correction to 
the pion-nucleon vertex contributes to photoproduction on the proton. At z/ = z/^m + 2 
there are tree-level contributions involving the magnetic moment of the nucleon. Together 
with the previous orders, this reproduces an old "low-energy theorem". However, still at 
^ = i^min + 2 there are large loop corrections where the photon interacts with a virtual 
pion which then rescatters on the nucleon before flying away to the detector. This quantum 
effect was missed in the old "low-energy theorem" and in most models. It is determined 
in terms of known parameters and thus constitutes a new, corrected low-energy theorem. 
At z/ = Umin + 3 other loops and new counterterms appear. Results for the cross-section at 
threshold up to z/ = z/^m + 3 from a fit to data above threshold constrained by resonance 



saturation of parameters ||55[ are satisfactory. To this order there is a prediction for the near 
threshold behavior of the 'yn — > ir^n reaction, which would be important to test. 

Compton scattering on the nucleon can also be calculated with this EFT. To z/ = z/mj„ + 2, 
the amplitude contains non-analytic terms and is completely predictive |^^. Expanded for 
small energies, it yields predictions for the polarizabilities, which are determined by the pion 
cloud. There is good agreement with data on the proton, and again it would be important 
to test the corresponding neutron predictions. 

We have seen that there are a number of tantalizing predictions from xPT concerning 
nucleon structure, particularly involving the neutron and neutral pions. These are hard to 
test because neutron targets and neutral pion beams do not come by easily... Nuclear targets 
and virtual pion effects can compensate for this, if we can formulate an EFT for nuclear 
systems that is consistent with xPT. We are thus naturally led to nuclear physics proper. 



IV. THE TWO-NUCLEON SYSTEM 

The two-nucleon system is where we hope to understand the basic modifications in the 
xPT power counting needed to tackle bound states. I will in turn consider the three regions 
of momenta that have been investigated: Q < Mnuc in Subsect. [IV A| , Q < M^j^ in Subsect. 



IV B| , and Q < Mqcd in Subsect. [IV C| , although this is not a historical order. 



A. Very low energies 

The lowest range of momenta serves to illustrate the method of EFT in the simplest 
non-perturbative context. Because only nucleons appear explicitly, we can go a long way 
following analytical arguments, while a good fraction of the physics of higher momenta is just 
a generalization of the same ideas to a more complicated situation, ^-matrix elements for 
NN scattering at momenta Q <C Mnuc can be obtained in zeroth order in a from an effective 
Lagrangian involving arbitrarily complicated operators of only A^ and its derivatives. The 
only symmetry constraints come from Lorentz transformations of small velocity, parity and 
time-reversal. Electromagnetic interactions can be included by adding all possible U{l)em- 
invariant terms. In the two-nucleon system it will give rise to photon exchange in addition 
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to dressings already present in one-nucleon systems. This will be discussed after we consider 
strong interactions in the a — > limit. 

The amplitude. The issue of establishing a power counting is largely independent 
of complications generated by spin and isospin, so for now consider [^ a system of two 



identical, stable spinless particles if) in their center-of-mass frame at energy E = k^/m — 
k'^/Am? + . . . much smaller than their mass m, their internal excitation energy, and the range 
of their interaction, which I will denote collectively by M. In addition to the kinetic terms 
in Eq. (p!9D, we have to consider the most general Lagrangian involving four nucleon fields, 
which can be written as 



o 

+i(C2-C^)V^Vv2(^V) + ---, (57) 

where the C2„'s are parameters of mass dimension —2(1 + n) that depend on the details of 
the dynamics of range ~ 1/M. The ". . ." stand for operators with more derivatives. 

Conservation of particle number reduces the T-matrix to a sum of bubble diagrams. The 
potential is given in momentum space by 

V{p, p') = Co + C2ip' +p") + 2C',p- p' + ..., (58) 

p {p ') being the relative momentum of the incoming (outgoing) particles. If Fourier- 
transformed to coordinate space, these bare interactions consist of delta functions and their 
derivatives. The bubbles require regularization and renormalization, which introduce a mass 
scale A through some cutoff function F{P / h?) [^. A bubble involving n > derivatives at 
the vertices gives 

hn{k) = j -^^l'^Go{kk)F{l'/A') 



(2vr) 
m 



27r2 



u=o ^ -''- 



(59) 



where On is a number that depends on the regularization scheme chosen, and R{x) is a 
regularization-scheme-dependent function that approaches a finite limit as x ^ 0. In general 
On and R{x) can be non- vanishing with either sign. The severe A-dependence comes from the 
region of high momenta that cannot be described correctly by the EFT; it can be removed 
by lumping these terms together with the unknown bare parameters into "renormalized" 
coefficients C2n . Only the latter are observable. The regulator-independent, non-analytic 



piece in Eq. (p9D, on the other hand, is characteristic of loops: it cannot be mocked up by 



W Although I refer to this as a cutoff, I am not assuming any particular regularization scheme. 
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re-shuffling contact interactions. From Eq. (|59D we see that we can estimate its effects by 
associating a factor mQ/An to each loop and a factor Q to each derivative at the vertices. 

The T-matrix is given by the diagrams in Fig. 0. For the first few terms one finds 
on-shell 
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Here dependence on the cut-off was ehminated by defining renormalized parameters C( 
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6*2 , and C2 from the bare parameters Co (A), C2(A), and C2(A): 



G 



(R) 



(Q 



(R)^ 



1 

MR) ~ 
^0 


1 


m 
27r2 


k + (- 

V 


_C2 


m 

47r2 


'i?(0) 
A 


2m2 



'2a 



+ 



Cn 4m2 



c;(^) 



+ ... 



C,+ 



(61) 



Other terms coming from R{k'^/A'^) are oc k'^ or higher powers, and cannot be separated 
from higher-order contact interactions that I have not written down exphcitly. 

It is important to reahze that Eq. (|60|) consists of two different expansions: a loop 
expansion and an expansion in the number of insertions of derivatives at the vertices or 
particle lines. The derivative expansion depends on the relative size of the coefflcients C2„ , 
for example C2 Q"^ /Cq . The loop expansion is governed by ttiCq Q/47i. The sizes of 
terms in both expansions thus depend on the absolute size of the 6*2^ s. If one makes the 
simplest naturalness assumption that there is only a single mass scale M and no fine-tuning 
in the underlying theory, we expect all parameters to scale with M according to their mass 
dimension. It is convenient, however, to factor in a 47i /m and write 6*2^ ~ 47i / mM'^'^^^ . 
In this case, the derivative expansion is in ~ (Q/M)'^, while the loop expansion is in ~ Q/M: 
the T-matrix in Eq. ( |60D is a simple expansion in powers of Q/M. However, in nuclear or 
molecular physics we are interested in shallow S-wave bound states ^ A shallow (real or 
virtual) bound state means that a parameter a of the underlying theory is close to its critical 
value Uc for which a bound state is exactly at threshold. That is, the underlying theory is 
so fine-tuned that it has two distinct scales: M and a lighter scale K = \a/ac — 1|M <^ M 



§§ The remaining case 
in tfie literature ISf 



-that of an unnaturally small scattering length — has also been considered 
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FIG. 2. First few terms of the two-particle amplitude T in a natural EFT. Two solid lines 
represent a Schrodinger propagator; a circle on a line represents a Q'^ relativistic correction; a heavy 
dot stands for a Q" contact interaction, and a dot within a circle for a Q^ contact interaction. 



induced by the fine-tuning. I will deal with an S'-wave shallow bound state by considering 
the class of theories where 



C. 



(R) 
2n 



Att 



m^M'ii)'^ 



72n, C. 



2n 



.{R) 



An 



^j^4-2n+l 



l2r. 



(62) 



with the 72n s dimensionless parameters of 0(1). The natural scenario is recovered when 
a is tuned out of Oc, and K becomes comparable to M. The loop and derivative expansions 
contained in Eq. ( |50D are in ~ Q/N and ~ Q^/KM. 

For Q <^ 'R, the theory is perturbative in QfR. As Q becomes comparable with K, all 
the terms in the loop expansion are of the same order {0{An/m'R)) while corrections from 
derivatives go in relative powers of K/M <^ 1 and can be accounted for systematically, at 
each order in K/M. The amplitude is (see Fig. ^ 

TUkp' ■ p) = TW(fc) + T«(A:) + Tj^fik) + ... + T^if^Kp' ■ p) + . . . 
= (Ti°)(A;))o - 2C'^^'^h^p'-p + 0iiATT/m^)i^/My). 

The S'-wave component can be rewritten, up to higher order terms, as 

iTos{k))o = [(Ti°)(A:))-i + T«(A:)(Ti°)(A:))-2 + TJ;f{k){TJ;f{k)r' + . . .]'' 
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(64) 
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FIG. 3. First three orders of the two-particle amphtude T in an EFT with a shahow bound 
state. Symbols as in Fig. ^. 



where C4 represents a certain combination of Q'^ contact interactions. By resumming the 
largest terms, we have produced a new expansion where the leading amplitude is 0(47r/m(K-|- 
Q)) and corrections go as Q'^/M{'R + Q) or similar combinations. The amplitude (Q) is 



nominally correct including 0((47r/'mK)(N^/M^)) only, but as we will see in a moment it 
is actually correct including 0((47r/mK)(K^/M^)). As we consider larger Q, Q ^ K, the 
leading term becomes OIAtt/itiQ) and corrections get relatively more important, growing in 
powers of Q/M. The new expansion fails at momenta ~ M, as expected. 



The amplitude (|63D and (|6^) is in the form of an effective range expansion in each partial 
wave 0,0, 
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In the S wave we obtain for the scattering length and the effective range 
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In the P wave, we find a scattering volume 

4" = ^. (68) 

Higher moments V2 can be obtained from the A;^" terms that arise in higher order in the 
Q/M expansion of the amplitude. For example, the S'-wave shape parameter is 

An S'-wave bound state can arise as a pole of the amplitude (B3) at k = in, 



«^ = — i^U + 2^(— i^l +0((K/Mf)|. (70) 



4vr ,Cr / 47r 



Cq > (< 0) implies k > (< 0) and represents a real (virtual) bound state. (The residue 
of i times the S'-matrix at this pole is indeed positive (negative) if k > (< 0).) The binding 
energy B2 to this accuracy is, of course, just the usual effective range relation among B2, a2 
and r2. 

We see that an /(> 0)-wave effective range parameter of mass dimension 5 is 
0(1/(2/ + 1)M^). The S wave is trickier. The scattering length 03 scales with K rather 
than M, 02 = 0(K~^). The S'-wave bound state is within the range of validity of the 
EFT, as K = 0(K) and the binding energy B2 = K^/m + Oi^X^/W^) < M. It is clear that 
the bound state can be treated in a systematic expansion in K/M. The two terms in the 
effective range rg (^) come at different orders. The main contribution originates in the 
contact interactions and scales with M, so that r^ = 0{M~^). The relativistic correction, 
on the other hand, is smaller by ~ M'R/m'^ ~ K/M; this justifies the neglect of relativistic 
corrections in many situations and the usefulness of a non-relativistic framework in shallow- 
bound-state problems. An underlying theory consisting of a non-relativistic potential can 
thus serve as a test for the scalings of 02 and r2 . Indeed, it has been explicitly shown |5^ 



that if the underlying theory consists of a non-relativistic potential of strength a and range 
1/M, then near critical binding effective range parameters behave precisely in way derived 
above. As for the shape parameter, it appears from Eq. (|69D that ^2 ~ 1/KM^, but 
such a scaling is seen neither in phenomenological analyses nor in models. This implies that 
74 = 7I/70 + 0()i/M): although each term in Eq. (HD comes at relative 0(Q^/M^K^), they 
are correlated and their sum is only of relative 0{Q^/M^'R). We conclude that the shape 
parameter appears at the same order as the first P-wave contribution when Q ^ N. Eq. 
(|6^) is therefore good up to 0(K^/M^), as advertised. This means that the first three orders 
of the expansion of the amplitude — depicted in Fig. H — are pure S wave and given solely 
by Cq and C2 ■ It is easy to verify that, in order to avoid M/K enhancements in higher 
effective range parameters, there must be such correlations in higher order coefficients as 
well: 

dJ) = d"^) (Cf ^ /Ci"^^ )" + smaller terms. (71) 
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FIG. 4. The dressed dimeron propagator: a bar is the bare dinieron propagator, while the 
bubbles represent two-particle propagation. 

We conclude that in the case of short-range interactions, the EFT approach, when applied 
consistently, is completely equivalent to the effective range expansion. The power counting 
described here and this equivalence were spelled out in a manifest scheme-independent way 



in Ref. ||5^, and in specific subtraction schemes in Refs. |60|j61[| . 

Dimeron. Although not necessary, it is possible to reorganize the EFT in order to 
account explicitly for the dominant correlations among the C^n s- Because they form a 
geometric series in C2 /C*o 5 ^^ can sum the dominant correlation terms as we have done 
above with the Cq contributions themselves. This procedure generates a four--^ interaction 
of the type Cg /(I — 2C2 k"^ /Cq ); the first two terms in Eq. (|6^) follow immediately. 
Now, this new resummation resembles the exchange of an s-channel particle. This justifies 



the suggestion [£2| that one can use a new EFT which involves, besides particles ip, also 
a "dimeron" field T with the quantum numbers of the shallow (real or virtual) two-body 
bound state. The leading terms in the most general Lagrangian consistent with the same 
symmetries as before are 

£ = aTt Udo + ^V' _ A + . . .) T - j={T^iPil, + h.c.) + ... (72) 



Here g and A are parameters and a is the sign of the kinetic term. The Lagrangian ([72| ) is 
equivalent to (^) if 

Thus after renormalization {g^^^f /2ti = 0{M/m^), A^^) = 0(MK/2m). ". . . " in Eq. (|72D 
stand for terms with more derivatives, which are suppressed by powers of K/M. Effects of 
non-derivative four-?/' contact term can be absorbed in g'^/A and terms with more derivatives. 
The coupling to two-particle states dresses the dimeron propagator. The dressed propa- 
gator contains bubbles as in Fig. ^ plus insertions of relativistic corrections. This amounts 
to a self-energy contribution proportional to the bubble integral. As we know, this integral 
is ultraviolet divergent and requires renormalization of the parameters of the Lagrangian 
([T^). Relativistic corrections can be accounted for as in the EFT without the dimeron. The 
dimeron propagator in its rest frame (p" = k'^/m — k'^/Am^ + . . .) has the form 

S{p^ ) = a^^^- ^jj^ . (74) 

The amplitude in this EFT is obtained directly from Eq. (|7^ — see Fig. |^. The result 
is identical to Eq. ( |64D to a given order. This is a more direct way to re-derive the effective 
range expansion, including the first few orders at once. The non-relativistic dimeron includes 
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FIG. 5. The two-particle amplitude T in the EFT with a dimeron field. 

the first two orders in K/M correctly: in leading order only the mass term in the dimeron 
kinetic terms needs to be retained, and observables depend on the combination (g^^-^Y / A^^-' ; 
in subleading order ^(^^^ and A*^^^ appear separately. The addition of relativistic corrections 
regains the first three orders. Corrections stemming from higher-derivative operators can be 
accounted for perturbatively, and will contribute to the ipip amplitude starting at relative 
0((N/M)^). The drawback is a mixing of different orders in the K/M expansion, with a 
number of smaller, higher-order terms included. 

The same real or virtual bound state as before exists near zero energy. The field T 
corresponds to a light (ghost or normal) state of mass A*^"^) ~ K relative to 2m. It serves 
to introduce through its bare propagator two bare poles at the k = ±y/mA. This is an 
undesired pole structure, but fortunately it changes character as the amplitude gets dressed, 
since the poles move to the imaginary axis. One ends up at the position of the shallow bound 
state, k = ing, k^ = k, given by Eq. (|70D, while the other pole signals the breakdown of the 
EFT at Kd = 0{M). The sign of the residue of —iS at each of the poles shows that the 
shallow pole corresponds to a regular (composite) particle, while the deep pole corresponds 
to a ghost. The field T can be used to describe the shallow bound state in processes other 
than i/ji/j scattering, if one keeps in mind that the residue of —iS at the pole is not 1, but 
{{g^^^y/g'^){a/{l + arn^g'^ /S'K^/^mB2)). Despite the resemblance, this is not an EFT with 
an "elementary" field for the bound state. At Q -C K, the EFT (with or without dimeron) 
can be matched onto a lower-energy EFT containing a heavy field for the bound state which 
does not couple to ipilj states, but the dimeron field shows structure at Q ^ ^^, which arises 
from if) loops. 

Schrodinger equation. Recall that an ordinary potential does not generate an 1/r 
behavior of the wave-function close to the origin because action of the Laplacian on such a 
singularity produces a delta function: V^(l/r) = — 47r(5(r). A local delta-function potential 
in turn is too singular to admit a solution of this type because then the potential side of the 
Schrodinger equation becomes V'^ ~ 6{r)/r, more singular than the kinetic side. This is 
the translation into Schrodinger language of the need for regularization and renormalization 
of the bubble diagrams. 

The effects of renormalization on the on-shell amplitude are two-fold. One is to replace 
the potential (|58| ) by 

!/(«) (k) = Cf^ + 2Cf )fc2 + . . . . (75) 

The other is to substitute the (free) Schrodinger propagator (|^) by 

Gf\r-k) = -^{l + 0{kr)). (76) 
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It can be shown |56| that this is equivalent to the usual Schrodinger propagation with 
an "effective renormalized potential" or pseudopotential. The renormalized Schrodinger 
equation is 

- -(V^ + k^)^{r) = -(Cf ^ + 2Cf ^fc2 + . . .)5(f)^(r*(r)). (77) 

fit CJ I 

Recalling the connections between EFT and effective range parameters, we see that this 
is nothing but a generalization of the well-known pseudopotential of Fermi and Breit 
||63H65|. Therefore, renormalization of an EFT with only short-range interactions leads to a 
Schrodinger equation with a pseudopotential determined by renormalized parameters, rather 
than a Schrodinger equation with a local, smeared delta-function-type potential determined 
by bare parameters. 

Further insight into the effect of renormalization can be gained by looking at the solution 
of Eq. (^). The effect of the -^r operator in the pseudopotential is to soften the singularity 
at the origin, replacing 5{r)/r by 5{r). Indeed, we can find the solution of Eq. (^) by 
expanding ^{f) = J2n=-oo^nr'^ at small r. The result is completely determined by A^i. 
Now, the same solution follows from 

1 Att 

(V^ + fc2)^(f) = —rdif^'^ir) (78) 

m m 

— which for any r 7^ is the free Schrodinger equation — with a peculiar boundary condition 
at the origin. 



47r /. ^ C 

r=0 



[K) 



mcf^r ^(ciV^---'- ^^^^ 



d 

ln(r\E'(f)) 
or 

We recognize the first term in Eq. ( [79D as the boundary condition of Bethe and Peierls, 
and Breit f^jQ' -^^^ renormalized EFT of only short-range interactions is equivalent to an 
energy expansion of the most general condition on the logarithmic derivative of the wave- 
function at the origin. The physics behind these Schrodinger formulations is clear if we look 
at the resulting radial wave-functions ripir) = u{r), which obey a free radial equation, (i) 
For P = -H? < 0, 

u{r) = C(_)e-'^^ (80) 

with C(_) a normalization factor and n given by Eq. (|TD|). (ii) For k"^ > 0, 

wW = C(+) sinfcr-- (o)"o),o,o , cosfcr , (81) 

y 1 — a2 r2 k^/^ + ■ ■ ■ J 

with C(+) a normalization factor and 02 and rg given in Eqs. (^) . By changing the 
unspecified dependence of C(+) on k and Ks,d we do not change the underlying pole structure 
of the ipip amplitude, but we do change the Jost functions. This is a reflection of the fact 
that the EFT cannot distinguish among short-range potentials which produce the same 
asymptotic behavior. Indeed, it is not difficult to see that with appropriate choices of C(+) 
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we can cover all situations considered in Ref. [|66|. The EFT shrinks to the origin the 
complicated wiggles that the "underlying" wave-function might have at small distances, 
thus replacing them with a smooth behavior that matches onto the tail of the underlying 
wave-function. The effective wave-function is all tail for the tail is all one can see from far 
away. 

Potential and regularization. I have so far avoided using specific regularization 
schemes in order to emphasize features of the renormalized theory, which is all that deter- 
mines observables anyway. Now a few words are spent on the related issues of the role of 
potentials and of regularization schemes. Some of these remarks have been presented in Ref. 
||56| , and shown to be supported by a careful numerical investigation [ |67[1 . 

It is clear that iteration of the potential is superfluous in any partial wave without fine- 
tuning: within its region of validity, the EFT is perturbative, and any bound states, if they 
exist, are beyond the range of the EFT. Iteration of the potential includes only part of the 
higher-order terms, and this incomplete set cannot a priori be taken seriously. It has been 
shown explicitly that the iteration of momentum-independent contact interactions does not 
reproduce the non-analytical behavior of a toy natural underlying theory ||6^ . On the other 
hand, since the error induced by iteration is small, it does not affect observables much. The 
more interesting cases are the ones with fine-tuning in the underlying theory. 

The appropriate formula for the amplitude, Eq. (0), was obtained by iterating the Q° 
part of the potential to all orders and by using perturbation theory for the corrections. The 
same result can be obtained more easily by first summing the bubbles with the full V{p',p) 
and then expanding in powers of the energy. For example, keeping track only of the Q^ 
terms in the S'-wave, 



(^"<*)'» = - (c^T2cW - ^'*))' • ^(*) = ^ + 2;^ ^'f*' - N + to?) 



27r2 
(82) 



This corresponds to the standard procedure of iterating the full potential, but it is now 
obvious that to remove the infinities we have to expand 



1 _ 1 0^2.2 

Co + 2C2k^ Co C? 



— - 2-|/c2 + . . . , (83) 



and consistently neglect Q'^ terms. We are then back to Eq. ( p^ truncated at 0()i./M). 
Terms of 0((5^) and higher can only be made regulator- independent by renormalizing higher- 
order parameters in the Lagrangian. Indeed, as we have seen before, O(Q^) terms in V{p,p') 
will contribute k^ terms to Eq. (0). 

Regulator-dependence is allowed insofar as it is small, that is, of the same magnitude 
as smaller terms that have been neglected. Cutoff dependence may arise from two sources. 
Some cutoff dependence comes from R{k'^/A'^), for it contains arbitrarily high powers of 
Q/A. Further cutoff dependence creeps in through bare parameters if we iterate the full 
potential to all orders but neglect to make the expansion (|8^) . For example, to 0((K/M)2), 
by solving the Schrodinger equation with a finite cut-off we induce a shape parameter 

^ mCl 7tA^\ ^ ' 2m2 ^ 7 8m^ ^ ' 
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We see that apart from a regularization-independent relativistic correction, cutoff depen- 
dence enters through both R and the bare parameter ratio CI/Cq. It is natural to ask 
what are the constraints on regularization procedures imposed by the requirement that 
both sources of cutoff dependence are small, as in more complicated situations it might be 
difficult to carry the step ( p3D explicitly. Note that the EFT theory with a dimeron field 
automatically performs the step (|83D , and therefore reduces differences among regularization 
schemes to the R they generate. 

In a cutoff regularization, the form of F{p / K^) is chosen explicitly. The error induced 
by R is removed order-by-order, the relative residual error being 0(K/A). Finite cut-offs 



have been strongly advocated in Refs. p5| , [70| . Keeping a cutoff A ~ M has one useful 
consequence. In this case, the large renormalized parameters C2„ can be generated by 
natural-size bare coefficients C2n = 0(47r/mM^""'"^). For example, 

Co(A) = -75^(1 + 7rK/27o^iA + . . .) (85) 

is of natural size even though Cq = 47r7o/mK is large. What a cutoff does is to give a 
scale for the bare parameters to be fine-tuned against. Before renormalization each loop 
contributes a large factor oc A but each Cq compensates with a small factor oc A^^, resulting 
a net ~ (1 + 7r(K + Q)/26iA + ...). The non-perturbative character of the theory is evident: 
it generates cancellations among the leading terms, and the relevant physics is contained 
in the terms that are suppressed by K/A or QfR. These remarks can be couched in a 
renormalization group analysis for the potential ||69|| . 

In this way the ratio of bare parameters can be made of natural size as well. The 
difference between a truncation of the rhs of Eq. ( p3D and the Ihs is then no bigger than other 
neglected small effects such as higher-order interactions in the Lagrangian: for example, the 
shape parameter (p^ ) induced by the cutoff is in this case 0(1/M^). A power expansion 
of the bare potential followed by its iteration to all orders is under these circumstances no 
worse than a perturbation treatment of the corrections to leading order. In this case the 
power counting (^5]) applies to the potential itself: since here the potential is a tree-level 
object involving four fermion fields 

u = d, (86) 

where d is the number of derivatives at the vertex. Truncating the potential at certain order 
accounts correctly for the first few terms in the expansion of the amplitude but includes also 
a subset of higher-order terms. We can actually invert the previous reasoning and use the 
freedom introduced by a finite cutoff to save some work. Given a regularization procedure, 
we can further fine-tune A to its "optimal" value that improves agreement with data. If we 
are dealing with a single channel, this is a consistent procedure: it simply means that 0{Q^) 
contributions have been included, but they do not need to be written explicitly because the 
corresponding bare parameter is zero. In general, however, only part of the higher order 
effects can be accounted for in this manner. As A — > oo the size of implicitly kept higher- 
order terms decreases. For example, all terms in Eq. (p^ ) — apart from the finite relativistic 
correction — then vanish. A consequence of this, noted in Ref. |^1|, is that the first two 



terms in V{p,p') are renormalizable in such non-perturbative context. 
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To any given order, equations relating bare and renormalized parameters will be trun- 
cated, and beyond leading order, be highly non-linear — see Eqs. (|6ll), for example. Inverting 
them, we can find the bare parameters in terms of the renormalized parameters (known from 
data) and the cutoff: C2n = C'2n(C'2n J ^)- This of course means that the bare parameters 
are not observables, depending on the scheme chosen. Given C2„ s, whether one finds real 
solutions C2n for any A depends on the actual values of the 6nS in Eq. (|59|) . In particular, 
to 0((K/M)^) the important term for large A is Cg O'^K^; the sign of Cg being that of 
■^ it is the sign of ^3 which is of concern. It was pointed out in Ref. |T^ that the most 



obvious cutoff schemes (such as a sharp momentum-space cutoff and a coordinate-space reg- 
ularization by square wells) imply 6'3 > 0. Imposing real bare parameters in such schemes, 
it is found that r2 < 2/A for large A, or in other words, that the limit A — i> 00 can only be 
taken if r2 < 0. This is a reformulation in field-theoretical language of Wigner's theorem 
| 73| : r2 cannot be positive if we take a sequence of ever shorter-range, (possibly non-local) 



hermitian potentials in a coordinate space Schrodinger equation. While such a constraint is 
unusual, it is doubtful that it is of much relevance to EFTs applied consistently to a certain 
order. First, this is obviously a regularization-scheme-dependent issue: it is easy to con- 



struct (smooth) cutoff schemes such that 6'3 < |74]. In the worst case scenario, the above 
constraint is a result of attempting to remove the cutoff in an inadequate regularization 
scheme. Second, this constraint only arises from a constraint imposed on unobservable bare 
parameters. No matter what process one is considering in the low-momentum regime, the 
bare parameters never appear, except in the right combinations with the cutoff that pro- 
duce renormalized parameters: an imaginary bare parameter is of no relevance. Indeed, the 
Schrodinger equation for the "renormalized potential" ( [77D violates the assumptions made 



in Ref. [^ about the nature of the potential. Finally, even if one insists on the unfortunate 
choice of regulator and the reality constraint on bare parameters, a positive r^ can still be 
achieved by keeping the cutoff finite and of the order of the mass scale of the underlying 
theory, A ~ M . 

An alternative regularization is dimensional regularization. F is implicit when we extend 
the dimension of space to a sufficiently small D — 1 such that the integral converges; the scale 
introduced is usually denoted /i. With minimal subtraction, 6'„ = and R{x) = 0. A non- 



minimal subtraction ("PDS") scheme was introduced in Ref. |60], in which a pole m. D = Z 
was subtracted. This particular scheme corresponds to 61 = 7r/2, 6'„>i = 0, and R{x) = 0. 
For n = this scheme reduces to minimal subtraction. It is not difficult to relate parameters 
in subtraction and cutoff schemes. If we call the cutoff parameter A = A(/i) + 7r/i/2^i, we 
can write 



2nK J 271-2 



_ n _ i,2(n+l) _ 

A Y: k^'r^n-^{K /i) + xA;2"(^fc + ^) + ^^R^e/K', /i/A) 
L i=0 ^ ^^ 



(87) 



where tt7„_j(A, /i) and R{k'^/A'^, fi/A) are related to the ^„'s and R{k'^ / A^). As we have seen, 
terms that grow with the cutoff can be cancelled by the counterterms C2n(A). If one defines 
new parameters C2n(/i) as A-independent functions of the C2n(A) and the ci7„_j(A,/i), we 
can express the renormalized quantities C2n as /x-independent functions of the C2n{f^) and 
fi. For example, in leading order 
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For /i, fc <^ A contributions from R are small; in the limit A ^ oo (with /i fixed) they vanish 
and the theory with cutoff can be rewritten in terms of the C2n(/^) and 

777 

l2n{k) = -—k^''{tk + fi), (89) 

which is equivalent to the theory with integrals subtracted at fc^ = — /i^ [|1^,^]. Other 



subtraction schemes in dimensional regularization have been discussed in Ref. [[75|1 . 

Dimensional regularization suffers from no errors generated by R. Minimal subtraction, 
first applied to this problem in Ref. ||7^, is special because it forces the bare parameters to 
coincide with the renormalized ones and be unnaturally large. Because no cutoff dependence 
is visible, if one works to Q^ terms it is not immediately apparent that l/(Co + 2C2/c^) has to 
be expanded in k'^. It was one of the results of Ref. [^ that a fitting of the ^5*0 NN phase 
shift based on the effective range expansion worked much better than one based on the form 
l/{Co + 2C2k'^). One can understand the failure of the l/(Co + 2C2/c^) fit by noticing that the 
induced spurious shape parameter (Rl) is here f 2 = 02 ^"2 /4 — 3r2 /Sm? + l/4a2 m'^, 



which is 0(1/KM^). Since there is no regulator to blame the fine-tuning on, the fine- 
tuning contaminates all other induced terms as well: they are all automatically large if 
02 is large. In dimensional regularization with minimal subtraction one cannot substitute 
the perturbative treatment of corrections by the iteration of the whole potential without 
reducing the range of applicability of the theory from M to the ratio of parameters in the 
power expansion of the contact interactions, which is a/KM. The perturbative treatment of 
corrections can still be carried out with dimensional regularization, or equivalently one can 
use the EFT with a dimeron field. This is true of minimal subtraction or other subtraction 
schemes. Of those, PDS and other schemes with a small subtraction scale are particularly 
convenient. For example, from Eqs. (|85|) and ( pS]) we see that natural-size bare coefficients 
translate into 

Coif^) = , ^l^, y (90) 

m{-fi + K/70) 

Co(K) = 0{4:7c/m'R) is a large running coupling. One can show that other C2„(/i) are also 
of the same size as their renormalized counterparts C2„ , and that they satisfy relations 
analogous to (0). We see that the power counting valid for the renormalized theory is 
manifest when the theory is formulated directly in terms of a small mass scale /i and running 
parameters C2„(/i) |]77| , |60| . 

Spin and isospin. Generalization to include spin S and isospin / is straightforward. 
In A^A^ scattering, S,I = 0,1. The Pauli principle requires / + S* + L be odd, where L is the 
orbital angular momentum. Channels are labeled by ^'^"'"^Lj, where J is the total angular 
momentum. There can be mixing between S = 1 channels with L = J — 1 and L = J + 1, 
generated by the tensor operator Sui^q) = 3q ■ aiq ■ (72/?^ — <?! ■ '?2- The corresponding 
amplitude can be expressed in terms of a mixing angle ej, 

.... _ .27r / e^*^-^-! cos 2ej - 1 ie^^^^-i+^-^+i) sin2ej \ . . 

^i^k))j.,,j+,--i—\ ^g.(5._,+5.+,)sin2£., e2^^^^+icos2£.;-l • ^^^> 
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One can use Fierz reordering to reduce the number of independent contact interactions 
to two Q° terms |10|, seven Q^ terms |12], and so on: in terms of g = p — fii and k = ^{p+ff), 



one possible choice is |T5| , p!6[] 

V{p,p') = d,'^ + C^PBi ■ a, + &^^q' + Cf P + (Cf )g2 + Ct^k')ar • a, 

^^^(5)^i_+^ ■ (g X A^) + Cf^q- aiq- 02 + cf ^A? • aiA? • 0^2 + • • • (92) 

The two Q" terms and two of the Q^ terms contribute to the two 5- wave channels, ^5*1 and 
^5*0. The remaining five Q^ terms contribute to ^5*1 — ^ D\ mixing and to the four P-wave 
channels, ^Pi, '^Pq, '^Pi, and ^P2. At higher powers of momenta, contributions to D and 
other waves appear. 

Since there are shallow bound states in both S'-wave channels, QCD must belong to the 



universality class of EFTs that obey the scaling (|62D with M = Mnuc- The only complication 
arises from the mixing generated by the Cg term. Analogously to the effective range 
expansion, we can expand 

1 sin2., = ^i^^= E u^^^''\kar) fc", (93) 

2 ^1 + k^a^^ n=o 

where the ^2 (^'^2 )'s are certain functions of the combination ka2 (but not of k 01 02 
separately). Models suggest that Wg (^^2 ) = ^(I/MImc"'')- This expansion is reproduced 
by the EFT provided 

with 7^ ^ dimensionless of 0(1). For example, uf^^^ oc Cf^ /{&q^ + C^P) ~ l/M^,. Such 
an argument can be generalized to higher orders. A local operator that connects a state with 
angular momentum / to another with angular momentum V involves at least / + /' derivatives 
and first appears at 0(((5/M„„c)''^'') if either / or V is zero, or at 0(((5/M««c)'^''^^) if both 
/ and /' are non-zero, relative to the lowest (0(47r/mK)) contribution (from / = /' = 0). 

Power counting then indicates that the low-energy phase shifts should be larger in the 5* 
waves, with the Q^ terms being treated non-perturbatively. The other waves, P, D, F, and 
so on, should be amenable to perturbation theory up to ever increasing energy, as we have 
seen that the suppression factor is ~ 1/(2/ + l)M„„c in an / wave. We also expect the waves 
to be roughly in decreasing order of magnitude as orbital angular momentum increases, as 
they receive contributions at increasingly higher orders. These conclusions are qualitatively 
substantiated by a PSA ||78| . 

Up to relative 0{Q/Mnuc) only 5* waves are present. In each S'-wave channel we have one 
non-derivative interaction at leading order and one two- derivative interaction at subleading 
order, which we can fit to the scattering length and to the effective range, respectively. 
From the singlet parameters a^p^°^ = —17.3 fm, a^p°^ = —23.75 fm, a^^°^ = —18.8 fm. 



r^^°^ = 2.85 fm, r^^"-* = 2.75 fm, and r^^''^ = 2.75 fm ||79| we find the averages 



^{Rrso) ^ ^3 . lo-3MeV"^ d'^K'So) ^ 45 . io-7MeV-^ (95) 
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FIG. 6. ^5*0 A'^A^ phase shift in radians as function of the center-of-mass momentum in the EFT 
to 0{Q/Mnuc) for several values of the cutoff /? (curves), compared to the Nijmegen phase shift 
analysis (circles). Figure courtesy of D. Phillips. 



From the triplet parameters a2 '^^' = 5.42 fm and r2 = 1.75 fm ||8^ we find 



A'Si) 



^(R)(35,) ^ 3 g . io-^MeV-2, Cf )('^^) = 2.1 ■ 10"«MeV-^ 



(96) 



Taking a channel-average K = 30 MeV and Mnuc = 140 MeV, these translate into reasonable 
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CSo) 



2.9, 
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eso) 



4.2. 



7r^" 



0.8, and 7^ 



i'Si) 



0.2. One can then predict the phase 



shifts as functions of energy. As an example, in Fig. ^ we see the resulting ^So phase shift as 
function of k for a number of cutoffs A = /3 |^ , compared to the values of the Nijmegen PSA 
| 78| . Cutoff dependence is significant only at momenta of 0(m,r) where this EFT should 
fail. The effective range expansion, which corresponds to the A ^ 00 limit, agrees well with 
the PSA up to such momenta. Similar conclusion holds for the ^Si phase shifts. In the 
case of the ^Si channel, we can alternatively use the deuteron binding energy to determine 
Cq , the fitting difference being of higher order. This is a more natural strategy for 

studying processes that involve the deuteron. 

At relative 0{{Q/Mnucf')-, mixing between ^Si and ^Di states appears thanks to C2 , 



and relativistic corrections need to be included in 5* waves 



(j(^){R) 



can be determined 

from the low-energy behavior of ei and then deuteron parameters can be predicted, or 
vice-versa. Using standard formulas in terms of integrals over wave-functions, deuteron 
properties found with C2 fitted to a deuteron asymptotic D / S ratio r] = 0.025 (as in 
the Argonne V18 potential [0) are given in Table |I| [Q. (Electromagnetic quantities refer 
to the contributions from lowest-order 'jNN couplings only, not to a consistent calculation 
which would include subleading one- and two-nucleon effects.) 



*The latter are numerically small, however, and can usually be neglected. 
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TABLE II. Results from EFT fits to subleading order for various cut-offs A and experimental 
values for the deuteron binding energy (B), magnetic dipole moment {fid)i electric quadrupole 
moment (Qe), asymptotic D/S ratio (77), and D-state probability (Pd). 



Deuteron 




A (MeV) 






quantities 


150 


216 


250 


Experiment 


B (MeV) 


1.799 


2.211 


2.389 


2.224579(9) 


ij-d im) 


0.868 


0.846 


0.828 


0.857406(1) 


Qe (fm2) 


0.231 


0.288 


0.305 


0.2859(3) 


rj 


0.025 (fit) 


0.025 (fit) 


0.025 (fit) 


0.0271(4) 


Pd (%) 


2.11 


5.89 


9.09 





The theory can be extended to higher orders in obvious fashion. 

Electromagnetic corrections. Because of the non-relativistic nature of heavy parti- 
cle propagation, electromagnetic interactions are dominated by Coulomb-photon exchange. 
One-Coulomb-photon exchange is 0(e^/Q^), while, using the power counting developed 
above, iterated once it is 0{me'^ /AttQ^). Therefore, Coulomb-photon exchange is (non-) 
perturbative for Q > am {Q < am). 

In the case of nucleons, amN — 7 MeV, which almost coincides with the light scale in 
the 5- wave channel where pp scatter, N*^ ^°'' ~ 8 MeV; that is, numerically, a ~ N^ '^"Y'^^at- 
For Q -C K*^ '^°\ Coulomb interactions have to be iterated to all orders, while we can treat 
strong interactions generated by the potential (p2| ) as perturbative insertions. Conversely, 
for Q ^"^^ '^"-', we can add electromagnetic interactions as perturbations. 

At Q ~ ^^*■ '^"\ both one- Coulomb-photon exchange and the momentum- independent 
contact interaction are of the same order and need to be iterated to all orders. Correc- 
tions can most simply be classified taking a ~ K*^ ^'^^ /Mnuc, and then recognizing that an 
extra numerical suppression affects electromagnetic corrections because Mnuc/^T^N -C 1- At 
(9(^( ^0) ^Mnuc), besides the two-derivative contact interaction, the most important electro- 
magnetic correction appears due to vacuum polarization, as it comes suppressed by a single 
power of a. (Actually, the more typical suppression of photon loops is a/ir, reducing fur- 
ther the impact of this effect.) At 0((K( ^'^'> / Mnuc)"^) , besides relativistic corrections to the 
S wave, transverse-photon exchange generates magnetic interactions. More exotic electro- 
magnetic interactions — such as two-photon exchange generated by the Thomson seagull — 



appear in higher orders. This ordering of electromagnetic interactions |^^ is in qualitative 
agreement with the findings of the Nijmegen PSA [^. 

The calculation of pp scattering including Coulomb interactions to all orders consists of 
dressing the bubble diagrams in Fig. ^ with Coulomb exchanges |8^j8^. This gives rise 
to integrals of the type (^) but with the free Schrodinger propagator ( ^5[ ) replaced by the 
Coulomb propagator 



Gocil;k) 



2Tram,iy/l 



rriN 



exp(arriAr//) — 1 P — k^ 



le 



(97) 



Dressing of the external legs gives rise to the usual Sommerfeld factor. In leading order, the 
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amplitude can still be summed as a geometric series in terms of a single bubble integral. 
In addition to the linear cutoff dependence present even in the absence of electromagnetic 
interactions, a new cutoff dependence am\i\n{K/am^)/TT appears. This requires a new 
leading-order, no-derivative four-nucleon interaction 

MR) _ 47ra .„„. 

to be added to Cg , which can be fitted to the difference a^p — a„„. In this EFT 

the splitting in scattering lengths cannot be calculated without matching to the underlying 
theory. Effective-range — but no vacuum polarization — effects were considered in Refs. |^6|| , 



where it was shown that to subleading order Vpp = Vrm- Effects of vacuum polarization and 
two-photon exchange are discussed in Refs. |87|j88|| . 

Moral. The modern EFT method applied to the problem of short-range forces is, after 
renormalization, completely equivalent to the ancient effective range expansion, and to the 
methods of pseudopotentials and boundary conditions. Why, then, bother with the EFT 
method? The gain is indeed marginal for the problem at hand, but there are advantages to 
a field-theoretical framework when other processes in the same energy scale are considered, 
and we want to treat them all consistently, free of off-shell ambiguities. (We will see an 
example in Sect. 0.) Yet the gain is potentially much more significant in more general 
cases where there are other, lighter degrees of freedom that generate longer-range forces. 
Symmetries are easily incorporated in Lagrangian field theory. The problem considered in 
this Section has very weak symmetry constraints. In more interesting problems, such as 
the NN system at energies comparable to the pion mass, the role of symmetry is more 
important, and the EFT method more convenient than the ancient techniques. We turn to 
this case next. 



B. Low energies 

As Q approaches m^r, it becomes increasingly difficult to account for pion exchange as 
a short-range effect. As we further increase momenta past Q ~ 171^^, we have to include an 
explicit pion field and build up all its interactions allowed by symmetries. The important 
new element is chiral symmetry, as discussed in Sect. PC| . Adding pions will generate all 
sorts of non-analytic contributions to the amplitude. For Q < m.„, pions need not but can be 
introduced explicitly, bringing non-analytic terms to k cot 6. However, these pionic terms can 
be expanded in a power series, and none of the terms with powers higher than the retained 
contact terms can be taken seriously, as neglected higher-order contact interactions could 
produce contributions of similar size. (Similar remarks hold for deuteron properties, as the 
deuteron is so dilute that it is even within the realm of the pionless EFT.) Nonetheless, one 
might hope to include at least part of the pion effects in this (uncontrolled) way. Indeed, low- 
energy properties seem to improve and cutoff dependence to decrease if one-pion exchange 
(OPE) is added to the contact potential at 0{Q/Mnuc) (0); which is then iterated to all 
orders [^. We want to go further and devise a rationale for a controlled expansion in 
the presence of pions, valid for momenta Q ~ Mnuc- The issue arises immediately of how 
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to estimate the size of pion contributions. One needs to find (i) tlie importance of pion 
exchange relative to short-range interactions; (ii) how iterated OPE compares to OPE. 

The OPE contribution to the A^A^ amphtude is, as we will see in more detail below, 
proportional to Q'^/M^^^{Q^ + M^^^). On the other hand, as we have just seen, lead- 
ing short-range effects are proportional to Att/Mqcd'^- Thus, as long as K < Mnn with 
Mnn ~ 47rM^^^/MQ(7D, OPE should be suppressed compared to the leading contact in- 
teraction by N/Mtvat. Using power counting arguments again, once-iterated pion exchange 
can give a contribution ~ (Q^/M^„^((5^ + M^^^))^(MQ(7£»Q/47r), and it also looks like a 
small correction to OPE as long as Q ^ AfjvAf. Of course, as in the lower-energy EFT, we 
expect all contributions to higher-/ partial waves to be perturbative at least up to momenta 
~ {2l + l)Mnuc- With such estimates, we conclude that at sufficiently small Q pions might be 
treated perturbatively. The suggestion that pions can be profitably treated as perturbative 
even in S waves and a power counting to support it were made in Ref. [^ in connection 



with a more complicated formahsm, but details have not been published. Ref. |^ gave the 
first clear implementation of these ideas. 

It is important to find out whether this can be done for momenta above 171^^. Only if Mnn 
^ m„ are perturbative pions an improvement over the effective range expansion. Moreover, 
the full power of chiral symmetry in constraining pion interactions relies on an ordering 
of chiral-symmetry-breaking effects in powers of m^. Including pions perturbatively brings 
terms that are a priori suppressed only by powers of rriT^/Mi^N- Chiral symmetry will be 
really useful only if this ratio is smaller than 1. The problem is that M^n is not obviously 
larger than Mnuc', whether or not it actually is depends on dimensionless numbers that cannot 
be accounted for with dimensional analysis alone. The issue of the range of validity of the 
EFT with perturbative pions can only be settled by explicit calculation of dimensionless 
factors and comparison with precise observed quantities, this being done to sufficiently high 
order that a significant number of pion effects is tested. 

Let us thus assume K < Mnn, Q < M^n, and ttIt, < M^n- A simple power counting 
can be done by taking Q ~ K ~ uIt, < M^n, and counting powers of the light scale Q. This 
is a direct extension of the power counting in the pionless EFT in Sect. [VA. In particular. 



the scahng of the contact operators is assumed to be the same as before with M = M^^, 
and thus their ordering is unchanged. Because of chiral symmetry, each insertion of a pion 
exchange brings a factor of Q/M^n- Electromagnetic interactions can be treated much as 
in the previous section. 

If MjsN is sufficiently high, we need to include the delta isobar as well. Beyond contact 
interactions considered previously in Eq. (p^), and the pion-nucleon interactions in the 
one-nucleon sector shown in Sect. Q, there are new interactions involving four fermion 
fields where at least one is a delta, and involving four fermion fields and pions. For example 
11,111, 

£(0) = -DrNhaN ■ {N^TSA + h.c.) + . . . , (99) 

£(^) = -DiN^N NH-g- DN - D2D ■ {NHaN x NHaN) + . . . , (100) 



TT^^ 



£(2) = mM 1 - 2D-'—\ {ci%N^NN^N + CglN^aN ■ N^ctN) + . . . (101) 
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FIG. 7. Insertions at 0{Q/Mi\in) in the perturbative-pion EFT, to be substituted for the circled 
dot in the second line of Fig. ^. A dashed line represents a pion, a heavy dot a derivative and a 
heavy dot with two circles two derivatives or powers of the pion mass. 



Here Dt, Di and 6*2^^ ^^^ undetermined constants of mass dimensions —2, —3 and —4 
respectively. ". . . " above stand for interactions involving more delta fields. So far investi- 
gations in this EFT have integrated out the delta isobar. As in the one-nucleon sector, this 
amounts to omitting the explicit-delta interactions and replacing the coefficients of nucle- 
onic interactions by new ones that scale with factors of l/5m. For example, Eq. ( p.00|) then 
involves new coefficients dj with 



d =D _^RlhA^ 
9 6171 



9 (5m 



(102) 



The amplitude. The leading-order NN amplitude coincides with that in the lower- 
momentum EFT, represented by the first line of Fig. ^ and by TJ^-* (the Cq and ik terms) 
in Eq. (Q): it is the iteration of chirally-symmetric, non-derivative contact interactions. At 
this order there are contributions only in the two S-wave channels. Subleading terms, of 
0{Q/M]yN) relative to leading, are constructed in a direct extension of subleading terms of 
the lower-momentum EFT, second line of Fig. ^ and Tj]) (the C2 term) in Eq. (|6^). The 
difference is that now besides two-derivative contact interactions we also insert one-pion 
exchange and a non-derivative contact interaction that breaks chiral symmetry explicitly. 
In other words, besides the C2 term in Eq. (|58D we need to consider one insertion of 



V = - 



f2g. 



\F^ 



' cri-qa2-q -^^ 

ti ■ *2^^-- — + m^[C. 

q^ + m^ 



(S) 
2qm 



+ C2qmP\ ■ ^2} 



(103) 



— see Fig. |^. Here we have neglected the energy transferred in OPE as it is a higher order 
effect. When Eq. ( p.03| ) is inserted within loops with TJ^-*, new divergences arise that are 
oc m^. These divergences can be absorbed in the C2gm's, giving rise to renormahzed C^qm^- 
Both C2 and C2qm, contact interactions also only contribute to S waves. The tensor operator 
from pion exchange, on the other hand, introduces mixing between ^5*1 and ^D\ waves. To 
this order, the ^5*1 — ^ D\ mixing parameter and the ^D\ phase shift are 'predicted in terms 
of pion parameters (and the short-ranged parameters that can be determined from the ^5*1 
phase shift). 



A calculation of the A^A^ system to this order was first carried out in Ref. pO] in the PDS 



subtraction scheme, and repeated with a coordinate space cutoff in Ref. [q^]. One finds in 
the ^5o channel [§g,||^ 



fccot(5('^o) 



,('5o) 



a^"^ 



+ 



A;2 + 5('^°)('^) + 0(A;VA,m2/A,l/4'^''^'A), (104) 
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This is in the form of an effective range expansion supplemented by a non-analytic "shape 
function" generated by pion exchange. Similar expressions hold for kcotS^ ^^\ and analytic 
formulas can also be derived for ei and 5^ ^^^ ||60|j89| -pT| . 

(i5o) ^(i5o) 



To this order, there are three unknown parameters in the ^5*0 channel (Cg ° , C< 



G 



('5o)^ 
2gm J 



and three (Q "^^ C. 



(^5i) ^(^5i) 



^2qm , 



in the '^S\ channel. Note, however, that the am- 



plitude to this order can be written as a function of Co + m?^C2qra for it is this combination 
that appears in the potential (p2D + ( p^03D . Only at higher orders will C^qm contribute new 
terms through pion loops. To subleading order, then, the separate dependence of observables 
on Co and C2qm is an artifact of perturbation theory. The difference between treating C2qm. 
perturbatively and non-perturbatively is an indication of the size of higher-order corrections. 

A fit to the Nijmegen S'-wave phase shifts [^] was done in Ref. [^ using the amplitude 
derived in the PDS scheme. As described in Sect. |IVA| , the theory can then be written 
in terms of running couplings such as C^qm. (/^) • The scattering lengths and effective ranges 
can be used to determine Co and C2 with a choice of C2qm- Choosing C^qmijn-K) = 0, the 
predicted ^5*0 phase shifts deviate significantly from the Nijmegen PSA at momenta of order 
m.^. The amplitude with three free parameters in each S'-wave gives a much better fit for 
momenta below 200 MeV; the predicted ^Si — ^ Di mixing parameter ei and the ^Di phase 
then come out reasonably well — see Fig. |^ |BD[ . 

Note that to this order other waves are determined as well. For example, P waves are 
given by OPE, and indeed it is well-known that the sizes of triplet scattering volumes is 
roughly reproduced by OPE [Q. Similarly, it is been part of nuclear folklore that high 
partial waves are mostly due to perturbative pion effects |^. The various perturbative 
contributions from pion exchanges to D and higher waves was examined in detail in Ref. 
9^ , up to sub-subleading order. The D wave cannot not be described well above laboratory 
energies of ~ 50 MeV without short-range interactions, but higher waves get progressively 
better described by perturbative pions, as expected. 

Schrodinger formulation. The EFT with perturbative pions can be recast in terms 
of the Schrodinger equation in coordinate space. In Ref. |^ it was shown that the above 
subleading amplitude can be obtained by a Born expansion on the Yukawa part of the pion 
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FIG. 8. ^5*0, ^Si, and ^Di NN phase shifts and ei mixing angle in degrees as functions 
of center-of-mass momentum in MeV: EFT in leading (dot-dashed) and subleading (dashed line) 
orders; and Nijmegen PSA (solid line). Figure courtesy of M. Savage. 



,Ui[r 



potential at distances r > R r^ 1/A. The wave-function can be written as u{r) = J2i 
where (i) the leading term is a solution of the free Schrodinger equation, uq = sin(A;r + 60), 
with k cot 6q determined by the contact interactions; and (ii) higher-order terms satisfy 
Ui{R) = 0. tii(r) can be expressed in terms of simple integrals involving the long-range part 
of the potential. 

Range of the theory and regularization issues. From the expressions for the 
amplitude shown above, it seems that pion effects are associated with the scale Mtvtv ~ 
iTiF^/g^niN — 300 MeV [Q. If so, the expansion parameter is ~ tti-^/Mnn — 0.5. The 
scattering length (|105|) is formally of the form ~ K + 0{Q'^ /Mj^j^i). Cq can be estimated 
in leading order using Eq. (|67|) and then used in Eqs. (|106|) and (|107|) to estimate the 
relative size of pion contributions. For example, OPE is responsible for about half of the 



effective range ( 106 ). Note, however, that this simple-minded procedure is not without its 
OPE contributions to Eq. 



caveats. OPE contributions to Eq. ( |105| ) are not small. They include a contribution which 
is non-analytic in the quark mass; this term clearly shows not only that the expansion is in 
m.„/Mj^]si even at /c = 0, but also that dimensionless numerical factors can be dangerous, 
since 5'^m^mAr/27rF^ ~ 0.9. Nevertheless, C2qm can be adjusted numerically to compensate, 
maintaining the lowest-order mostly unchanged. The reasonable success observed in the low- 
partial-wave phase shifts suggests that dimensionless numbers do indeed conspire to increase 
MjsN- A number of recent papers have attempted to provide a more precise estimate oiMj^^q. 
The issue of a practical test of convergence has been tied to specific regularization schemes 
employed in calculations. 

In Ref. [P^l the ^Sq fit has been reexamined with both PDS and a subtraction at p^ = 
—fi'^. Significant differences were found in the phases for Q ^ fi depending on whether C2 
and C2qm were included perturbatively or non-perturbatively. This led the author of Ref. 
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to conclude that pions can be included perturbatively only for external momenta well 
below the pion mass. In Ref. 



95 



a similar conclusion was reached through a Lepage error 
analysis. In Ref. |9^ it is argued that this strong sensitivity to /i arises from the fact that 
the subleading amplitude in PDS is not /i independent, so that the parameters that satisfy 
renormalization group equations cannot give a good fit to phase shifts. Pion exchange gives 
contributions that do not vanish in the chiral limit, and introduces subleading but potentially 
significant dependence on /x. It was shown that this fi dependence can be ameliorated by 
splitting Cq into two pieces, Cq = Co[o] +Co[2], where Co[o] is treated non-perturbatively while 
Cq[2] ~ rnlC2gm ^ C*o[o] is treated as a first-order correction. The amplitude in subleading 
order is then a function of Co[2] + rn1C2qm- An improved subtraction scheme ( "OS" ) designed 
to minimize the /i dependence of the perturbative expansion was also introduced, and the 
subleading amplitude computed. Fits to ^5*0 and ^Si phase shifts for momenta between 7 
and 100 MeV — fits that are weighted towards low momenta — are then reasonable, at the 
same time producing good values for the scattering lengths. The freedom of fitting C2qm 
can enforce smallness of the full subleading contribution. The authors of Ref. ||96| conclude 
that there are no obstructions to using perturbative pions for external momenta larger than 
the pion mass. 

Refs. I^J^ point out that the shape functions are completely determined by pion ex- 
change. For momenta below the pion mass, the shape functions can be expanded in powers 
of fc^, and effective range parameters other than scattering lengths and effective ranges are 
likewise determined. For example, in either S'-wave channel ^2 ~ '^/f^^^'^-^NN- Such 
relations can be viewed as a one-pion-saturation model [^ of the parameters of the very- 
low-energy EFT of Sect. |1V A| . One can see that the perturbative-pion EFT provides no 
explanation for the fine-tuning in the scattering lengths, but does provide higher-order effec- 
tive range parameters which scale approximately as 1/M^~^ as expected from dimensional 
analysis. Note that these parameters blow up in the chiral limit as a consequence of the long 
range of pion interactions. If the perturbative-pion EFT can be shown to hold for Q ^ m.,^, 
then these predictions can be elevated to bona fide low-energy theorems. On the other hand, 
if the perturbative-pion EFT holds only for Q ^ m,r; contact parameters showing up at the 
same order render the pion predictions no better than an order-of-magnitude estimate. The 
shape function and higher effective range parameters look like good testing grounds of the 
convergence of the expansion as — contrary to C2qm in ^2 — there are no parameters that 
can be fine-tuned to balance OPE effects. 

In Ref. |^D[ the shape functions were expanded in a series k"^. The mixing parameter ei 
was expanded as well. 



£i 



Y.a^'^k 



2i+l 



(108) 



although this is expected to breakdown at very low momenta ~ K; in view of Eq. ( P3| ) it is 
no surprise that gi ~ l/K^*~^m7rM7VAr. The result of a comparison between perturbative pion 
predictions and values of effective range (|66D and Si ( |108|) parameters extracted from the 
Nijmegen PSA |j7^ is in Table [Q. Because of the dramatic failure of these predictions 
the authors of Ref. 0] suggest that the scale of the physics left out from this perturbative- 
pion EFT is not much larger than irij^. They further suggest that successes of this EFT 
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TABLE III. Parameters of analytic expansions of kcotS in the ^5*0 and "^Si channels and of 
the ^Si —^ Di mixing angle ei predicted by the EFT with perturbative pions compared to values 
obtained from the Nijmegen PSA. 
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are coining from the expansion in 'R/M^n, with the failure of the expansion in m.,r/MNN 
disguised by counterterms. They attribute the good fits of Fig. |^ to the fact that any 
theory with free parameters consistent with the effective range expansion can reproduce the 
rapid rise at low k and the subsequent slow decrease of the S-wave phase shifts. However, 
using the energy-independent Nijmegen analysis, the authors of Ref. [Q have argued that 
the uncertainty in f 2 is too large to make a definite test of the OPE predictions. They 
claim that uncertainties in the ^2 s from the energy-dependent Nijmegen analysis are 
unrealistically small, and question whether low-energy data alone can be used to determine 
these parameters. Ref. |^T| has reexamined these issues directly in terms of the shape 



function at finite energy. It is shown that the shape function in the ^Si channel is well- 
determined by the PSA: at the deuteron pole, the shape function is close to zero to a good 
precision resulting from the accurately known values of B2, 02, and r2; at positive energies 
the finite-energy phase shifts are also need to extract S. It is found that OPE saturation 



fails again way beyond error bars. Because S is small compared to k cot 6, Ref. [^ points 
out that one has to include relativistic and isospin-breaking effects in the extraction of S. 
This is a little unsettling, however, as these are higher-order effects that were not included 
in the EFT estimate. 

Although work done so far has gone some distance in assessing the validity of perturbative 
pions, this remains one of the most important open questions for the EFT approach in 
nuclear systems. Work along this line requires the evaluation of higher-order effects. At 
sub-subleading order we encounter new pion exchanges: both non-static (or radiative) OPE 
and once-iterated OPE. Radiative pion diagrams have very recently been studied in Ref. ^7_ 



It is suggested that the size of these graphs can be estimated counting powers of y/m^m^, yet 
mysterious cancellations have been found that render the sum of leading non- vanishing non- 
static pion diagrams smaller than expected. The role of deltas has not yet been examined. 
Particularly in diagrams with non-static pions, integrating deltas out should significantly 
limit the range of the EFT. So far there is also no explicit calculation of the dimensionless 
factors that are involved in once-iterated OPE, and of the sizes of sub-subleading contact 
interactions. Work is in progress in this direction PSI and a better determination of the 



range of validity of this EFT should be available soon. As we will see in Sect. [VA| , this is 
pivotal in assessing the usefulness of this EFT for nuclei other than the deuteron. In Sects. 
[VA| and |VIA| it will be evident that this EFT allows a much simpler treatment of nuclear 



systems than afforded by traditional methods. 



45 



Finally, note — in complete analogy to the extensive discussion in Sect. |IVA| — that 
if pions can be treated perturbatively then little error is made in iterating them to all 
orders. There might even be something to gain if one is lucky that short-range effects are 
small because of unknown dimensionless factors. The procedure of including pions in a 
potential that is then iterated to all orders had actually been carried out before any of the 
developments reviewed in the present section. We turn to this next. 

C. Moderate energies 

We want to extend the EFT approach to momenta Q ^ M^n- How this is to be done 
depends critically on the value of M^vat. At this scale pions become non-perturbative. If 
MjVAT is as large as Mqcd the task is daunting as there is presently no rationale for a 
controlled expansion in the presence of massive degrees of freedom such as the rho. On the 
other hand, if Mjvtv -C Mqcd, we might hope to improve the expansion of the previous 
section by a controlled resummation of selected terms that go as Q /Mnn- If M^j^ < 
300 MeV, the lack of known particle thresholds there suggests that the resummation could 
involve only pions. For Mattv ^ 300 MeV delta isobars probably need to be included as well. 

In fact, naive dimensional analysis alone suggests that the leading effect of pion exchange 
should be iterated to all orders for Q ~ Mnuc- This is because numerically Mqcd ~ 47rM„„c, 
and thus M^n ~ Mnuc- The new element is that OPE is ~ 1/M^„^ ~ An/MQCDMnuc- 
Because OPE has short-range components, it is natural to assume that they set the scale 
for K, K ~ Mnuc "C Mqcd- Part of the unnaturalness of K is thus explained. (It is 
still mysterious why K < Mnuc, some cancellation between short and long-range effects is 
necessary to enforce this.) This suggests a power counting might be meaningful in which 
we take Q ~ M]^is[ ~ K ~ Mnuc ~ Mqcd/4:71. Now, arguments analogous to those in Sect. 



IV A| get complicated to the point that they can no longer be carried out analytically in 
detail. Nevertheless they are still valid: a loop containing a reducible intermediate state 
still contributes MqcdQ/'^t^, so adding no-derivative contact interactions (~ An /Mqcd"^) 
or static OPE (~ 1/M^^^) to a graph amounts to a power of Q /)^ ~ QjMnuc- In order 
to describe Q ~ M^uc we have to sum all graphs with the leading potential, no-derivative 
contact interactions plus static OPE. All other contributions would come as corrections. This 



power counting is due to Weinberg [|T0|,|T3|, but many of its elements had been anticipated 



by Friar PIJTOOl. 



The question has been raised, whether such a power counting can be carried out consis- 
tently. The problem is that once OPE is iterated to all orders, renormalization of the theory 
becomes quite different from the perturbative-pion EFT. Perturbative arguments have been 
given ||60|j96[] to support the claim that Weinberg's power counting is not consistent. Ref. 



|60| points out that a diagram with OPE between two Co interactions requires the C^qr, 



counterterm, as we have seen in the previous section. Ref. |^6[ finds that thrice-iterated 
OPE and diagrams with twice-iterated OPE plus one Cq interaction require a C2 countert- 
erm. This implies that alleged higher-order contact interactions might get contaminated 
by powers of Mqcd I Mnuc in the numerator, to the point that they are required to appear 
in leading order. If this is a feature of an infinite number of contact interactions, there is 
no ordering and this EFT is doomed. However, it has been known in the context of the 
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Schrodinger equation that perturbative arguments of this type are not in general rehable 
for singular potentials, as is the case for OPE. The perturbative expansion might have a 
cut starting at g\/ F^ = 0; insistence on a g\/F^ expansion would then reflect itself on 
different orders offering correlated contributions to counterterms, each bringing powers of 
MgcD/Mnuc, yet resulting in a much better behaved sum. This conjecture ultimately has 
to be decided on the basis of actual calculations. We will assume it to be correct, fit phase 
shifts and then examine the cutoff dependence and convergence of the expansion. If the 
conjecture is correct, then we should find that this EFT has a range of validity greater than 
Mj^ii^. Although less familiar to those raised with perturbative calculations, this procedure 
is quite analogous and no less justified than the methods being used to assess the range 
of validity of the perturbative-pion EFT. We are going to see that there is evidence that 
renormalization is being performed correctly, as cutoff dependence turns out to be weak. 

The power counting is most easily done for the potential V. An arbitrary contribution 
to the two-nucleon potential can be represented by a two-nucleon irreducible Feynman di- 
agram, that is, with A = 2 continuous nucleon lines and at least one pion or delta isobar 
in intermediate states. Effects from the scale Mqcd are lumped in the vertices, and by 
restricting ourselves to the potential we exclude the scale Q"^ /Mqcd from energy denomi- 
nators. Since the only explicit scale is Q, the xPT counting (M) applies directly if we put 



A = 2 mM 



iy = 2L + J2V,A,. (109) 

i 

In this way the benefits of chiral symmetry are extended to nuclear physics: since both L 
and Aj are bounded from below {L > 0, Ai > 0) , u > Umm = 0. Note that this counting is 
a direct generalization of the counting (^) of the pionless EFT, where the potential has no 
loops and a single vertex of index A = d. Isospin violation will be considered shortly. 

We assume implicitly that once Mn^c {F^,, ttIt,, possibly 6m, and, neglecting fine-tuning, 
K) have been accounted for, all other dimensionful quantities are given by Mqcd- Since we 
also assume that no enhancements appear from the iteration of the potential, the scaling 
of contact interactions is determined by pion exchange in the potential: as in the case of a 
single nucleon, we expect additional Q^ to be accompanied by Mq^jjj, so that 

with the 72n s dimensionless factors of 0(1). 

Thanks to infrared enhancement, reducible diagrams have to be treated differently than 
the potential. The leading-order potential is iterated to all orders. Corrections to the leading 
potential do not need to be iterated to all orders. Yet, as we have seen in Sect. [IV A| , they 



can be iterated with small error as long as one uses a regularization with a cutoff mass 
A ~ Mqcd] for the two-nucleon system, this is equivalent to solving the Schrodinger equation 
with the potential V. A potential up to a certain order ensures that the amplitude is correct 
to the same order. Of course, even if our assumptions about the scaling of interactions in 
this EFT are incorrect, results from the iteration of the potential with a sufficient number of 
interactions automatically include all results of perturbative pions discussed in the previous 
section. 
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Traditionally, potential models have been plagued by problems of principle, such as 
the form of meson-nucleon interactions (for example pseudoscalar vs. pseudovector pion 
coupling), renormalization issues, absence of a small expansion parameter, etc. Because 
the EFT potential includes explicitly the exchange of only pions, all these problems can be 
resolved. As we have seen in Sect. [11 C| , because pions are pseudo-Goldstone bosons, the form 
of their interactions with nucleons is determined by the pattern of chiral symmetry and its 
breaking. (For example, if the EFT were formulated relativistically, pseudovector coupling 
would be preferred, as insistence on pseudoscalar form would demand new seagull vertices 
to ensure "pair suppression" P3). Renormalization can be performed as all interactions 



consistent with symmetries have been included. And the power counting ( |109| ) for the EFT 
potential implies that diagrams with an increasing number of loops L — and, in particular, 
with increasing number of exchanged pions — should be progressively less important. Clearly, 
the EFT potential can be thought of as a low-energy approximation to standard potential 
models, although this can only be taken in an average sense. The scope of an EFT potential 
for systems with at least two mass scales has been examined in toy models. Ref. |7^ 



considers the case of a short-range potential supplemented by a Coulomb potential. It 
exemplifies through an efficient error analysis how the approximation of contact interactions 
for the short-range potential improves the low-energy description systematically. The small 
role played by the cutoff A can also be shown explicitly in this case. Similar conclusions hold 
in toy models where two nucleons interact via a light scalar with Yukawa coupling and a 
short-range potential [p5| , |101|| . It has been shown that, even in the presence of a fine-tuning 



that produces a shallow S'-wave bound state, the contact interactions obey Eq. (|110|) , with 



Mqcd representing the heavy mass scale. (Ref. |^ further emphasizes the advantages of 
a fit to modified effective range parameters.) The EFT potential thus provides a QCD 
justification for some of the phenomenological successes of potential models. 

Potential and amplitude. The Lagrangian for this EFT is the same considered before 
in the context of perturbative pions. A calculation of all contributions to the two-nucleon po- 



tential up to 1/ = I'min + 'i was Carried out in Refs. HT^JTSUTql using time-ordered perturbation 
theory. Some diagrams are shown in Fig. ^. 

In leading order, u = Vmini the potential is simply static OPE and momentum- 
independent contact terms, that is, the first term in Eq. (|103|) and the first two in Eq. 



(p^) [|T^. This is obviously a very crude approximation to the NN potential. It does al- 
low us to see how energies of 0(10 MeV) arise from QCD: the average of the leading-order 
potential in coordinate space is roughly (V^V) ~ —g\m^/AT^F^ ~ —10 MeV. However, it is 
known that the nuclear force has other sizable components, like a spin-orbit force, a strong 
short-range repulsion and an intermediate range attraction. These are all generated in the 
next orders: v = Vmin + 1 corrections vanish due to parity and time-reversal invariance, but 
V = ^min + 2 corrections are several. First, there are short-range corrections; they come 
from one-loop pion dressing of the lowest-order contact interactions, and from four-nucleon 
contact interactions with two derivatives or two powers of the pion mass, displayed in Eqs. 
(0) and (|103|) . It is easy to show that the result of loop diagrams amount to a simple shift of 



the contact parameters. In Refs. |T2|,|T5|,|l6[ the Cq's were redefined in order to reabsorb both 
loop contributions and the chiral-symmetry-breaking C2gm's, as to this order they cannot 
be separately determined from A^A^ data. Second, there are corrections to OPE; these come 
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FIG. 9. Some time-ordered diagrams contributing to the two-nucleon potential Va^a^ in the 
EFT. (Double) solid lines represent nucleons (and/or deltas), dashed lines pious, a heavy dot an 
interaction in C^^\ a dot within a circle an interaction in C^^\ and a dot within two circles an 



interaction in C^"^'. First line corresponds to 



second and third lines to v 



+ 2, 



fourth line to v 



+ 3, and 



denote v > I'min + 4. All orderings with at least one 



pion or delta in intermediate states are included. Not shown are diagrams contributing only to 
renormalization of parameters. 



from vertex dressing and from recoil upon pion emission: 



V, 



(2) 



2g, 



tree 



F2 
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(Ti-q(T2-q 



Aiq^ + A^k"" - 2g, 



E-T^M^' + r 



q,2 _j_ ^2 



;iii) 



Again, here loop corrections have been absorbed in mass and coupling constant renormaliza- 
tion. Also, chiral-symmetry-breaking corrections to the tcNN coupling proportional to m^ 
can be lumped in gA_, which then satisfies the Goldberger-Treiman relation without discrep- 
ancy. Third, there are two-pion exchange (TPE) diagrams built out of lowest-order nNN 

and ttA^A interactions of Eq. (^7]). Denoting u± = \j{qi: ly + Am^^ these loop diagrams 

. g-2 _ p ^^^2 _ f2)2^ ^^ . ^^^ f)^^ . ^^^ f)) 



can be expressed schematically as 

(2) _ {9A,hA) ,. . . . f dH 
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At V 
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;ii2) 



Vr, 



, + 3 a few more TPE diagrams appear, which involve the ttttNN seagull vertices 



from Eq. (^8]); again schematically. 



V, 
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loop 



{9% h\)B, 

Fi 
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dH K (g2 _ r2)^ (^-2 _ f2)2^ ^^ . ^^ ^ f)^^ . (^-. ^ f)) 



(2vr)^ 



(a;±,ci;± 5m) 



■ (113) 
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To this order there are some relativistic corrections as well, but they are numerically small 
and have been neglected. 

This potential has all the spin-isospin structure of phenomenological models, but its 
profile is determined by explicit degrees of freedom, symmetries, and power counting. The 
power counting suggests a hierarchy of short-range effects: 5* waves should depend strongly 
on the short-range parameters Cq] contact interactions affect P-wave phase shifts only in 
subleading order, so their effect should be smaller and approximately linear; D waves are 
affected by contact interactions only via mixing, while higher waves should be essentially 
determined by pion exchange. Chiral symmetry is particular influential in the TPE piece. 
The latter includes a particular form of terms previously considered by Brueckener and 
Watson [|102|| , Sugawara and von Hippel [|103|| , and Sugawara and Okubo ||104|| , plus a few 



new terms. Those terms involving the Pj's and the deltas provide the only form of correlated 
TPE to this order, as graphs where pions interact in flight appear only in next order and 
should thus be relatively small. The sum of the Bi term and the corresponding delta term 
(which give the C3 of Eq. (|54D, related to the nucleon axial polarizability) , is particularly 
important in providing an isoscalar central force. Not surprisingly, in the chiral limit these 
potentials behave at large separations as van der Waals forces. 

Regularization and renormalization are necessary not only for the loops in the potential 
but also for the loops generated by the solution of the Schrodinger equation. For numerical 
convenience a smooth cutoff of the Gaussian type was used, and calculations performed with 
the cutoff parameter A taking values 500, 780 and 1000 MeV. For each cutoff value a set 
of (bare) parameters was found that fits phase shifts below 100 MeV laboratory energies. 
There are nine independent parameters stemming from contact interactions, although for 
the fit to A^A^ phase shifts we have also varied nine other, redundant parameters. OPE 
and TPE diagrams are completely determined by parameters accessible in pion-nucleon 
reactions, but because most were not known at the time they were also searched in the fit. 
A sample of the results for the lower, more important partial waves is presented in Fig. 
|10|, together with phases from the Nijmegen PSA |7^ — cf. Fig. §. Quality of the fits is 



typical of other waves. Waves higher than F were found to be mostly described well by 
pion exchange alone, as expected. Deuteron quantities are shown in Table ^ — cf. Table 
|l[ (Electromagnetic quantities refer to the contributions from lowest-order 'jNN couplings 
only, not to a consistent calculation which would include sub-leading one- and two-nucleon 
effects.) The predicted S'-wave scattering lengths (not used to constrain the fit) were found 
to be 02 ° — —15.0 fm and a2 — 5.46 fm. Important central attraction comes from 
the Sj's and deltas, and indeed the central potential does resemble that from models that 
include a and u meson exchange explicitly. Values for the parameters are listed in Ref. ||16|| . 
Reasonable values were found for quantities known at the time, for example qa = 1-33 (in 
agreement with the Goldberger-Treiman relation) and h^ = 2.03 (smaller but not too far off 
the large- AC value). However, the values for the -Bj's came out different than those found 
later in vriV scattering: using Eq. (||), ci = 2.2 GeV"\ C3 = -0.7 GeV"\ and C4 = 2.7 
GeV~^, while from the leading delta contribution alone (because B4 appears here only in 
higher order) C2 ~ 1.8 GeV^^. These values are to be compared to those in Table |; one 
concludes that B^ is too big and of wrong sign, while B2 is too small and Bi too big, and/or 
the delta contributions are too small. For a cutoff of A = 780 MeV, the coefficients C2n of 
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FIG. 10. ^5*0, ^5*1, anii ^D\ NN phase shifts and ei mixing angle in degrees as functions of 
the laboratory energy in MeV: EFT up to z^ = 3 for cutoffs of 500 (dotted), 780 (dashed), and 1000 
MeV (soHd line); and Nijmegen PSA (squares). 



the contact interactions were found to scale approximately as in Eq. ( |11(]| ) with Mqcd ~ 500 



MeV for Mnuc = F^. 

A number of papers have subsequently examined different aspects of this EFT. In Ref. 
||105|| it was pointed out that the parameters Ai in Eq. ( |111| ) needed not be included explic- 
itly in the fit, as one combination is fixed as a numerically small relativistic correction ||106|| , 
while the Ai contribution can be absorbed in the C2 of Eq. (|92D and in the tiNN coupling. 



Note that the above EFT potential is energy dependent. As pointed out in Sect. piD| , equiv- 
alent potentials can be obtained through unitary transformations. An energy-independent 
potential is more convenient in many situations, and the corresponding version was derived 
in Refs. [ 107 , 105 , 106|| . The effect of the unitary transformation is essentially to remove the 
third term in Eq. ( |111D and at the same time to modify the TPE potential. It corresponds 
in what box diagrams are concerned to a change from the Brueckener and Watson ||102| 
to the TMO ||108|| potential. The latter has a further advantage: it displays explicitly a 



51 



TABLE IV. Results from EFT fits at i^ = 3 for various cut-offs A and experimental values for 
the deuteron binding energy (B), magnetic dipole moment (fid), electric quadrupole moment (Qe), 
asymptotic D/S ratio (ry), and D-state probability (Pd)- 



Deuteron 




A (MeV) 






quantities 


500 


780 


1000 


Experiment 


B (MeV) 


2.15 


2.24 


2.18 


2.224579(9) 


ij-d im) 


0.863 


0.863 


0.866 


0.857406(1) 


Qe (fm2) 


0.246 


0.249 


0.237 


0.2859(3) 


rj 


0.0229 


0.0244 


0.0230 


0.0271(4) 


Pd (%) 


2.98 


2.86 


2.40 





cancellation of the isoscalar spin-independent central and of the isovector spin- dependent 
interactions, which with a energy-dependent choice comes only through a combination of 
terms in the potential and its iterations. The components of the EFT TPE potential were 
studied in more detail in Ref. ||9^. In particular, it is shown explicitly that (i) relativistic 
corrections are mostly small (although in the weaker components they can make relatively 
large contributions); (ii) both isoscalar central and spin-orbit potentials are numerically sim- 
ilar to a and u exchange in models, if the i?j's take values given by ttN scattering; (iii) the 
OPE isovector tensor potential is reduced by the TPE contribution; (iv) a subset — which 
is invariant under pion field redefinition — of the u = Vmin + 4, two-loop TPE diagrams is 
small and repulsive. (Other modern studies of chiral aspects of TPE include Ref. ||109|| .) 

The fair agreement of this first calculation and data up to center-of-mass momenta 
Q ~ 300 MeV suggests that this may eventually become an alternative to other, more 
model-dependent approaches to the two-nucleon problem. In particular, it hints that a a 
exchange might be unnecessary. The relative cutoff insensibility also might be indication 
that renormalization is being performed correctly. The large values of the cutoff and the 
reasonable value of the contact terms are auspicious for the range of validity of the the- 
ory. Note that to this order this EFT includes all interactions incorporated to date in the 
perturbative-pion EFT and then some. More detailed analyses of the effects of different 
components of this potential on the fit have appeared since. 



70 the same EFT 



In Ref. 

potential (with the same type of cutoff) as in Refs. |l^JT5|JT6|| but without TPE (and energy- 
dependence in OPE) was refitted to the Nijmegen PSA |^ in three channels, ^5*0, ^Pi, and 
^D2-, below center-of-mass energy of 100 MeV. An error analysis shows explicitly how the 
second-order corrections improve over the leading-order fits, particularly if in a global fit 
weights are assigned to data. It is confirmed that pion effects contribute little to the ^5*0 
phase shift in this energy range (as expected from fine-tuning), but play a more important 
role in higher- wave phases (as expected by power counting). The error analysis hints at a 
breakdown of the expansion at M ~ 300 MeV, and shows that cutoff effects are minimal 
when A ~ M. It is suggested that the breakdown at M ~ 300 MeV is related to the ne- 
glect of TPE. Similar conclusions were reached in Ref. |^ . More recently, (the tail of the 
energy-independent version of) the EFT TPE potential in the limit of a heavy delta was 
substituted for (the tail of) the one-boson exchange in a Nijmegen phase-shift reanalysis 
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of pp data below 350 MeV ||110||. A drop in x^ was observed, ci was found to be poorly- 



determined by the fit; fixing it from the sigma-term determined by the modern vrA^ PSA of 



Ref. [ |111|| , the best-fit values are 

ci = -0.76(7) GeV"\ Cg = -5.08(28) GeV-\ C4 = 4.70(70) GeV'^ (114) 

The agreement with vrA^ scattering values in Table | is remarkable and confirms unequivocally 
the validity of chiral TPE. Taken together, these new results suggest that a more extensive 
search of parameters might improve the quality of the EFT potential fit. 

Isospin violation. We can order isospin-breaking effects in the potential in the same 
way done in Sect. |ITl| ||r5| , [T8[| . At Q ~ Mnuc 3> K, photon exchanges are indeed perturbative, 



as discussed in Sect. |1V A| . Here I will concentrate on effects beyond those considered there, 
that is, on those stemming from the quark masses, from indirect electromagnetic effects, 
and from simultaneous pion-photon exchange. I follow the standard nomenclature and refer 
to an isospin-symmetric potential as "class I" , to a potential that breaks charge dependence 
but not charge symmetry — defined as a rotation of vr around the 2-axis in isospin space — 
as "class 11" , to one that breaks charge symmetry but vanishes in the np system as "class 
III" , and to one that breaks charge symmetry but is present in the np system as "class IV" . 
No isospin-violating effects appear at leading order, z/ = Vmim so the leading potential is 
class I. The first effect appears at z/ = Vmm + 1 in the form of a class II isospin violation from 
pion mass splitting (Am^ oc aMqfjj;)) in OPE and from Coulomb photon exchange. One 
order down, v = Umin + 2, a class III force appears mainly from the quark mass difference 
through breaking in the ttNN coupling {f3i = 0{6m'^/Mq^^)) in OPE, from contact terms 
{.1s,a = C){em1/MQQ^)), and from the nucleon mass difference (Am at = 0{eml./MQCD))- 
Therefore, to order u = Vmin + 2 the isospin-violating nuclear potential is a two-nucleon 
potential of the form 

V = Vu[itiUt2h - tl ■ t2] + lill[(tl)3 + (^2)3], (115) 

where 

2gApiq-aiq-a2 , ^ ^. , . 

Vm = —=;- ^5- - (7, + 7^ai ■ ^2). (117) 

P:; q^ + m^ 

Finally, class IV forces appear only at order u = z/mm + 3. One concludes that the symmetries 
of QCD naturally suggest a hierarchy of classes in the nuclear potential ||T5|,|T8|: 



(^m+i)/(^m) ~ 0{Q/Mqcd), (118) 

where (Vm) denotes the average contribution of the leading class M potential. This qual- 
itatively explains, for example, the observed isospin structure of the Coulomb-corrected 
scattering lengths, a„p ~ 4 x ((a„„ + app)/2 - anp) ~ 4^ x (opp - a„„). 
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FIG. 11. Electromagnetic loops contributing to isospin violation in the one-pion-range nuclear 
force in first non-vanishing order. Solid lines are nucleons, dashed lines are pions, and wavy lines 
are photons. Other orderings are also included, but iteration of static OPE and Coulomb exchange 
is excluded from (e). 



One can use the above formalism to do consistent and systematic calculations of isospin 
violation. The isospin- violating potential of range ~ l/^Tj-Tr up to z/ = Umin + 3 was com- 
puted in Refs. | J112| , p.liJ| |. It consists of diagrams where one pion flies between the two 
nucleons, constructed out of tree-level diagrams with interactions of index up to 3, plus 16 
(not counting trivial particle permutations) non-vanishing one-loop diagrams, comprising 
pion and nucleon self-energy corrections, pion-nucleon vertex corrections, and simultaneous 
pion-photon exchange 



-see Fig. |TI 



Most of these diagrams are ultraviolet divergent; di- 
mensional regularization was used and l/(-D — 4) terms cancelled against bare coefficients. 
All diagrams with photons are also individually infrared divergent; a photon mass is intro- 
duced and removed at the end, the sum being infrared finite. The result is invariant under 
both gauge transformation and pion-field redefinition: 



Vny{q) 



7rF2 



ti ■ to 



(il)3(t 



2)3) 



(Ji ■qa2-q 

g2 



2n2 



(1 - gymj) 
[1 + q^/mDq^/ml 



ln(l + q^/ml 



(119) 



Its isospin structure allows only charged-pion exchange and therefore affects only np scatter- 
ing. This 7r7 potential has been incorporated in a Nijmegen phase shift reanalysis of np data 
below 350 MeV |p.l2|| . We can use the values for the iiNN coupling constants determined by 
the analysis to find that their isospin breaking is consistent with zero, with an uncertainty 
comparable to our expectation from dimensional analysis and from it — t] — rj' mixing ||114 . 
Similarly, 7^ might be viewed as originating in p — u mixing, while pseudovector-meson 
exchange (in particular close- lying doublets such as cti — /i) exchange contribute to the 7^ 
spin-spin force |114||115| . 



V. FEW-NUCLEON SYSTEMS 

Few-body systems provide non-trivial testground for ideas developed in two-body dy- 
namics. As we have seen, the deuteron is sufficiently dilute that an EFT for Q ~ ^^ can 
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be used to describe it well. Such an EFT is essentially equivalent to an effective range 
expansion. The latter has been successfully applied to reactions involving the deuteron and 
photons, but difficulties have been encountered in the three-body system. In Sect. |V A| we 
will see how the EFT can be extended to the three-body system. In this case some fairly 
unusual renormalization takes place, which makes the EFT method much more challeng- 
ing and thus interesting than in its more standard perturbative applications. The EFT for 
higher momenta where pions are treated perturbatively shares the same leading order as 
the pionless EFT, differing only in that pion exchange appears explicitly together with a 
contact correction in subleading order. The important question to be answered is whether 
few-nucleon bound states can be dealt with in the approximation where leading-order inter- 
actions are taken as short ranged. In Sect. |VB| still higher momenta are considered, and 



the connection to potential model approaches discussed. We can use this EFT to organize 
the various few-body-force structures. 

A. Very low and low^ energies 

Symmetries do not forbid three-nucleon contact interactions. In order to extend the 
EFT to the three-nucleon system we need to find the size of three-nucleon forces relative to 
two-nucleon forces. As in Sect. [IV A| , I omit at first spin and isospin variables. I also neglect 



electromagnetic interactions. Requiring invariance under small- velocity Lorentz, parity, and 
time-reversal transformations, at sufficiently low momenta, the Lagrangian involving six 
fields ijj is 

/: = -:^(V^W + ..., (120) 

o 

where ". . . " stand for terms with more derivatives, which are suppressed at low momenta. 
Note that an S'-wave contact three-body interaction is allowed for bosons. The same is true 
in the case of three nucleons in a spin S = 1/2 state, but the Pauli principle forbids an S'-wave 
interaction for three nucleons in a S = 3/2 state. The first possible three-body interaction 
in the latter channel contains an extra Q^, and should be relatively less important. 

A comparison between the tree-level connected three-body diagrams allows us to guess 
the relative size of two- and three-body interactions. We saw in Sect. [IV A| that operators 



involving four nucleon fields contain the small scale K. At momenta Q ~ ^^, two sequential 
two-body Co interactions contribute 0((47r/r?2N)^m/K^). If we assume that the coefficients 
D2n of the terms with 2n derivatives in Eq. ([120[) are of natural size, D2n ~ (47r)^/mM^+^"', 
then we expect three-body-force effects to be of relative ©((K/M)^"*"^"). But does the small 
scale K not contaminate operators involving six nucleon fields? This question is intimately 
related to the renormalization group flow of three-body interactions with the mass scale 
introduced in the regularization procedure. This flow, in turn, depends on the behavior 
of the sum of two-body contact contributions to three-body amplitude as function of the 
renormalization scale, or equivalently, as function of the ultraviolet cutoff A. 



It is convenient to rewrite this theory by introducing the dimeron field T of Sect. [IV A 
It is straightforward to show — for example, by a Gaussian path integration — that the 
Lagrangian ( [120[) is equivalent to Eq. (^) supplemented with 
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FIG. 12. The amplitude T for particle/bound-state scattering as a sum of dressed pinball and 
three- body- force diagrams (first and second hues) and as an integral equation (third line). 



C = hT^Ti)^ + ... 



121) 



In leading order, since the kinetic term of the dimeron can be neglected, observables will 
depend on the parameter appearing explicitly in Eq. ( |121| ) only through the combination 
-3 V/A2 = Do. 

Let us consider elastic particle/bound-state scattering at Q ~ ^^ | |116| - |11^| . (Three- 
body bound states manifest themselves as poles at negative energy. Both the inelastic 
channel and three-particle scattering involve the same type of diagrams considered here 
and can be obtained by a straightforward extension of the following arguments.) Diagrams 
contributing to the amplitude are depicted in Fig. 
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We have already seen that at these 
momenta all terms in the full dimeron propagator of Fig. ^ are equally important. It is 
easy to see that an L-loop diagram involving only two-body interactions ( "pinball" diagram) 
is Oi^img^ IQ"^) x {mg^Q/4:Tr{A + mg'^Q/ Att))^). Pinball diagrams are thus (superficially) 
ultraviolet finite, and are all of the same order for mg'^Q/AnA = Q/K ~ 1. The three- 
body amplitude is then the solution of an integral equation. This equation, including the 



three-body force, is also depicted in Fig. 12 



I choose the following kinematics: the incoming particle and bound state are on-shell 
with four-momenta {k'^/2m, —k ) and {k'^/4m — B2, k), respectively. The outgoing particle 
and bound state are off-shell with four-momenta {k'^/2m—e, —p) and {k'^/4m—B2+e, p); the 
on-shell point has e = k'^/2m — p'^/2m and p = k. The total energy is E = 3k'^/Am — B2. At 
very low momenta, we can limit ourselves to the more important S-wave. Using the boosted 
dimeron propagator (73), after performing the integration over the fourth component of the 
loop momentum, setting e = k"^ /2m — p^ /2m, and projecting onto the S-wave, we find that 
the amplitude a{k,p) normalized so that a(0, 0) = —03 is the particle/bound-state scattering 
length satisfies ||118|| 



a{k,p) = M{k,p]k) + 



2\ r^ 



vr JO 



dq M{q,p; k)- 



A;2 



le 



-a{k,q), 



(122) 



with the kernel 
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M{q,p-k) 




1 , fq'^ + qp + p^ 
— m 
pq 



mE 



h 



q2 _ qp _^ p2 _ jyi^j^ 



mg^ 



(123) 



Here A = 1 for bosons. The same equation is valid in the case of fermions, with different 
values of A. For three nucleons in a spin S = 3/2 state, this equation holds with A = —1/2 
and h = 0. For three nucleons in a spin S = 1/2 state a pair of coupled integral equations 
applies, with properties similar to the bosonic case. Eqs. (|122|) and (|123|) reduce to the 
expressions found in Refs. |11£, 116| , 117|| when h = 0. Note that the perturbative series 



shown in Fig. [12] corresponds to a perturbative solution of the integral equation for small 
A. The solution of Eq. ( p.22| ) is complex even below the threshold for the breakup of the 
two-particle bound state due to the ie prescription. To facilitate the discussion one uses 



the function K{k,p) that satisfies the same Eq. ( p,22| ) as a{k,p) but with the ie substituted 
by a principal value prescription. K{k,p) is real below the breakup threshold. a{k, k) and, 
consequently, the scattering matrix can be obtained from K{k,p) through 



a{k, k) 



K{k,k) 



l-ikK{k,k) 



(124) 



In order to understand the renormalization group flow due to two-body forces, let us first 
take h = 0. When p ^ l/a2 (but k ~ 1/0,2), the inhomogeneous term is small (0(l/pa2)), 
the main contribution to the integral comes from the region g ~ p, and the amplitude 
satisfies the approximate equation 



K{k,p) 



4A 



dq 



In 



q'^ + pq + p^ 



\^iT Jo q Kq^-pq + p"^ 



K{k,q). 



(125) 



Now, in the limit A — i> 00 there is no scale left. Scale invariance suggests solutions of the 
form K{k,p) = p"^, which exist only if s satisfies 



8A sin^ 

D_ 

\/3s cos f 



0. 



(126) 



If K{k,p) is a solution, K{k,p1/p) is also a solution for arbitrary p^,. Because of this special 
conformal symmetry, the solutions of Eq. ( |126| ) come in pairs. 

For A < Ac = 3v^/47r ~ 0.4135, Eq. (jT^) has only real roots. For A = -1/2, there are 
solutions s = ±2, ±2.17, ... In this case an acceptable solution of Eq. ( p,25| ) decreases in the 
ultraviolet, 



K{k,p:$> l/as) =Cp' 



(127) 



For finite cutoff, the solution should still have this form as long as p <^ A. The overall 
constant C(A) cannot be fixed from the homogeneous asymptotic equation, but is determined 
by matching the asymptotic solution onto the solution at p ~ l/a2 of the full Eq. ( |122| ). 
Because of the fast ultraviolet convergence, the full solution to Eq. (|122|) is expected to 
be insensitive to the choice of regulator, so that a well-defined limA^ooC'(A) can be found. 
There is thus no evidence for a contamination of three-body forces by the small scale K 
in the S = 3/2 channel of 5'-wave nucleon-deuteron (Nd) scattering. Up to and including 
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FIG. 13. k cot5 in fni ^ for S'-wave Nd scattering in the S = 3/2 channel as a function of k'^ in 
fm-2 in the EFT to 0(1) (dashed hne) and 0{{i</Mnucf) (soHd hne). Circles are from the PSA 



in Ref. [|120| and the triangle is from Ref. |121[ . 



relative 0{{Q/M)^), only two-body forces need to be included. A calculation using Eq. ( |72D 
up to and including relative 0{{Q/M)^) was carried out in Refs. ||116| , |117|| . As expected, 
cutoff dependence was found only for p ~ A. The result for the phase shifts for energies up 
to the break-up point is shown in Fig. |T3|: the result of the leading-order (0(1)) calculation 
| |119| | is given by the dashed line and the 0((K/M)^) result [ p. 17] ] by the solid line. These are 
compared to data points at finite energy from a PSA of Nd scattering data 



1201, and to 



the much more precise (nearly) zero-energy point ||121 
to be of 0(K/M„„c)3). At zero p-— -• — ^--i -^^/^^ 

6.35 ±0.02 fm 



We expect errors in our calculation 
ITTBI 



6.33 ± 0.10 fm 



compared to 



121| |. This is an extremely accurate low- 



the experimental value 03 
energy theorem! Our results seem to deviate from a simple effective range type expansion 
only around the pole at ~ 0.05 fm~^ . (A pole in k cot6 corresponds to a zero in the 
scattering matrix, which does not carry special meaning.) This pole does not appear in 
potential model calculations, and presumably will be (re) moved by higher-order terms that 
we have not yet included. The calculation of higher-order corrections involves the knowledge 
of further counterterms like the ones giving rise to S — D mixing. These parameters can 
be determined either by fitting A^A^ data or by matching with another effective theory — 
involving explicit pions — valid up to higher energies. If more precise experimental data 
— particularly at zero-energy — appear, we would be facing a unique situation where high- 
precision calculations in strong-interaction physics can be carried out systematically and 
tested. 

For A = 1, on the other hand, there are purely imaginary solutions: s = ±zso, where 
So ~ 1.0064. The solution of Eq. ([125D for p < A is 

P 



K{k,p ^ 1/02) = C cos SqIu 



A. 



(128) 



where, on dimensional grounds, A^(A) oc A. This solution has two constants C(A) and 
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A^(A), and implies that the solution of Eq. (|122|) is not unique in the limit A ^ oo 
||122|| . The undetermined phase A^ arises from special conformal symmetry, as there is 
no way here to select a preferred oscillatory solution. Solutions of Eq. ( |122| ) for A = 1 



and finite A (but h = 0) were obtained in Ref. ||118|| . We have verified that in the region 
1/(22 ^ p ^ A the solutions indeed have the form ( |128| ), with A^ approximately linear in the 
cutoff. The amplitude C was found to be approximately proportional to cos(soln(A^a2)). 
This dependence generates some points where the amplitude is cutoff independent, but 
apart from them, the amplitude does not show signs of converging as A — > cxd. Since the 
solution for small p has to match onto the large-p solution, the cutoff dependence leaks 
into the low-momentum region. Small differences in the asymptotic phase lead to large 
differences in, for example, the particle/bound-state scattering length. This leakage of high- 
momentum behavior into the low-momentum physics is indication that we are not performing 
renormalization consistently with our expansion. Note that if one were to truncate the series 
of diagrams in Fig. ^ at some finite number of loops one would miss the asymptotic behavior 
of K{k,p) that generates this cutoff dependence. This is because Sq (and its expansion in 



powers of A) vanish in a neighborhood of A = 0. The truncation of the series in Fig. |T^ is 
equivalent to perturbation theory in A, and cannot produce a non- vanishing sq. We are here 
facing truly non-perturbative aspects of renormalization, which resemble a phase transition. 
Faced with regulator dependence, one needs to modify the leading-order calculation, 
that is, to change the bookkeeping of the higher-energy behavior of the theory through the 
addition of at least one new interaction. In an EFT, one adds local counterterms. Here 
the only alternative is to introduce a three-body force, but the present case is complicated 
by the fact that the cutoff dependence of the amplitude is non-analytic around p = 0. 
This, however, does not mean that the renormalization program in this low-energy EFT is 
doomed: a three-body force term of sufficient strength contributes not only at tree level, but 
also in loops dressed by any number of two-particle interactions. It is convenient to rewrite 
the three-body force as 

h{k) = ^^. (129) 

H{A) has to be at least big enough to give a non- negligible contribution in the p r^ k r^ 
A region. This means that the dimensionless quantity H{A) has to be at least of 0(1). 
Assuming minimal strength, the three-body force has the feature that its contribution to the 
inhomogeneous term is small compared to the contribution from the two-body interaction, 
as it is at most p/A of the latter. When p ^ l/a2 (but k ~ l/ct2), the amplitude satisfies 
a new approximate equation, with a three-body interaction added to the rhs of Eq. ([1251) . 



The solution K{k,p ~ A) has a complicated form. But, because in the range l/a2 <C p <C A 
the three-body force term is suppressed by p/A compared to the logarithm, the behavior 
( |128| ) is still correct in the intermediate region. Now the phase is function of both A and 
H{A), so if H{A) is chosen appropriately, we can make A^ cutoff independent. Matching 
with the p ~ l/a2 solution should then determine the scattering length 03 = 03 (A^) and the 
low-energy dependence of the amplitude. For p ~ 1/02, Eq. (|122|) becomes 

3 K{k,p) 1, fp'^+pk + k^ -ME\ 



4 ("x + J3P1 _ ^e) P^ \P^ -pk + k^-MEj 
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FIG. 14. k cot(5 in units of l/a2 for 5- wave particle/bound-state scattering in the S" = case 
as a function of k in units of l/a2 in the EFT to leading order, for A^ = 19. 5a^ and different 
cutoffs (A = 42.6, 100.0, 230.0, 959.0 x a^^). 
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K{k,q), 



(130) 



where /x is an arbitrary scale such that /i <^ A, and we have dropped terms which are smaller 
by powers of p/A or yu/A. One can use this equation to obtain an approximate expression 
for H = H{A). By integrating out modes between A' > A and A and demanding that the 



low-energy amplitude be unchanged, we find |118 



H{A) 



sin(goln(£) - arctg(^)) 
sin(soln(£) + arctg(^))' 



(131) 



12| can absorb the high-momentum modes of the diagram that 



The last diagram in Fig. 
precedes it. 

The full solution of Eq. (|122D now depends, apart from two-body parameters, on the 
new, physical parameter A^ (but not on the regulator). This equation was solved numerically 
with a non-vanishing H{A) for several cutoffs. The three-body force that is necessary to 



keep a^ unchanged was found in good accord with the analytical approximation ( |131| ). In 
Fig. 14 we see the results for the corresponding K{k, k)~^ = kcotS, where 6 is the S-wave 
phase shift for particle/bound-state scattering, for different cutoffs in the case 03 = 1.56a2. 
The effective range, for example, is predicted as r^ = 0.57a2. 

We can extend the preceding analysis to the three-body bound-state problem. Since 
the inhomogeneous term did not play an important role in the above ultraviolet arguments, 
the latter hold for any bound state with binding energy B3 comparable with l/ma^. Such 
bound states are shallow, having a size ~ \j ^JmBl, ~ 02, and thus should be within range 
of the EFT. In principle, all bound states with size larger than ~ r2 should be amenable to 
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FIG. 15. Three-body binding energies B-^ in units of l/ma2 as functions of the cutoff A in units 
of l/a2 in the EFT to leading order in the S = case, for A^ = 19. 5a^ . 



this EFT description. Their properties will be determined in first approximation by only Co 
and A^, while more precise information can be obtained in an expansion in N/M. In Fig. 
T5| binding energies are shown for a range of cutoffs, with the three-body force adjusted to 
give a fixed scattering length 03 = 1.56a2. (For particular values of 03 and A, these results 
reduce to those in Ref. |[I23|| .) As we can see, for this value of A^ there exists a shallow 



bound state at B^ ~ 2/ma\ whose binding energy is independent of the cutoff. This bound 
state has size ~ 0.7a2 and can thus be studied with the EFT. As we increase the cutoff, 
deeper bound states appear at A„ = Aoexp(n7r/so) with n an integer and Aq — 10, so that 
for A.„_i < A < An there are n + 1 bound states. Similar results can be obtained for other 
values of 03. By repeating this procedure we can in fact determine 03 = 03 (A^) (or vice- 



versa] 



and S3 = S3 (A, 



or vice-versa] 



Eliminating A^, a universal curve S3 = S3 (03) can 



be obtained. In the case of the shallowest bound state, this is seen in Fig. 16. This curve is 



known in the three- nucleon case as the Phillips line ||124|| . It has been derived before in the 
context of models for the two-nucleon potential that differ in their high-momentum behavior. 
Varying among two-particle potential models one could expect to fill up the S3 x 03 plane. 
In the EFT, the high-momentum behavior of the two-particle potential is butchered, which 
causes no trouble in describing two-particle scattering at low-energies, but in the peculiar 
way described here, does require a three-body force. Varying among two-particle potential 
models is thus equivalent to varying the one parameter A^ of the three-body force. This 
spans a single curve S3 = S3 (03) in the S3 x 03 plane. Our argument suggests that the 
Phillips line is a generic consequence of renormalization group invariance. 

Corrections to this calculation come from operators with more derivatives, and they are 
discussed in Ref. ||118|| . The first correction comes from C2 in Eq. ([57|). It can be shown that 
the large A dependence it introduces can be absorbed in h itself, so that a finite, calculable 



contribution remains. Its effects are being calculated ||125|| . It is reasonable to assume then 
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FIG. 16. Binding energy i?3 of the shallowest bound state in units of l/ma^ as function of the 
particle/bound-state scattering length 03 in units of 02 in the EFT to leading order in the 5 = 



case. 



that more-derivative three-body forces do not get further enhanced by inverse powers of K, 
so that 



D 



2(n+l) 



D 



■(R) 



2n 



M2 



(132) 



In this case effects of three-body-force corrections could be accounted for perturbatively. 

The above results hold for ^He atoms | 118 |. In fact, it has recently been established that 
the two-'^He bound state ("dimer") is very shallow, with an average size (r) = 62±10 A ||126|| , 
more than an order of magnitude larger than the range of the interatomic potential. The low- 
energy two-^He system should then be describable by contact two-body forces. In leading 
order, the measured size translates into a scattering length 02 = 124 A, which determines the 
strength of the contact interaction. Unfortunately, although the three-^He ( "trimer" ) has 
been observed ||127|| , there seems to be no low-energy information on its properties nor on 
^He/dimer scattering. At least one three-body datum is needed to determine A^. Until such 
datum becomes available, we can only illustrate the method by using a phenomenological 
^He-'^He potential as a model of a microscopic theory. We select a potential [p,28|| that is 
consistent with the recent measurement of the dimer binding energy. It gives for the two- 
body system 02 = 124.7 A and r2 ~ 7.4 A. Three-body calculations are much more difficult to 
perform with such a phenomenological potential. An estimate for the ^He/dimer scattering 
length is 03 = 195 A. Ground and excited bound states have been reported; estimates for 
the shallowest bound state place it in the range -B3 = 1.04 — 1.7 mK, while a deeper state lies 
around S3 = 0.082 — 0.1173 K. There exists a prediction for the low-energy S'-wave phase 
shift, albeit for a different potential, but r^ could not be determined. Using such model 
we can estimate the range of validity of the EFT in momentum to be ~ l/r2 ~ 0.14 A~^, 
and the leading order to give an accuracy of ~ T2/a2 — 0.06, or about 10%, at momenta 
Q ~ l/a2. Using this potential's 03/02 = 1.56 we can determine A^, and predict both the 
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energy dependence in ^He/dimer scattering and the trimer binding energy. In fact, in Fig. 
p!4| I have used exactly this value of 03, so it represents the lowest-order prediction for k cot S 
for atom/dimer scattering, from which we can extract an effective range r^, = 71 A. From 
Fig. |l^ we predict a bound state at B3 = 1.2 mK, which is certainly within the EFT. The 
next-to-shallowest bound state is small enough to be at the border of EFT applicability. For 
a sufficiently large cutoff, we find B3 = 0.057 K, but in the best case scenario corrections 
from higher orders should be very important. 

Because of the similarity of the integral equations, our arguments should apply to systems 
of three fermions with internal quantum numbers as well. Preliminary results for the three- 
nucleon system in the J = 1/2 channel confirm this expectation and produce good Nd 
scattering and triton results ||125|| . The results of this section hold in the EFT where the 
pion has been integrated out, and can be straightforwardly generalized to the EFT with 



perturbative pions. Work in the latter is now in progress for the J = 3/2 channel ||129|] . We 



also plan to extend these calculations to the four- and more-nucleon systems ||125|| . We have 



seen in this section how we need to introduce a three-body force in leading order to prevent 
the collapse of the three-body system that would result from purely attractive two-body 
contact interactions. Is this repulsion enough to allow a prediction of other few-body bound 
states? 



B. Moderate energies 

Increasing the range of the EFT, we affect the ultraviolet behavior of the leading-order 
two-body amplitude. As we hit Q ~ Mjy^, contributions from pions become comparable to 
the leading contact interaction. The results of the previous subsection do not necessarily 
apply, as the two-body interaction is no longer purely attractive, but might display short- 
range repulsion besides long-range attraction. Here I continue in the vein of Sect. [IV C 



assuming that no enhancements appear in contact interactions. The power counting of that 
section can then be extended straightforwardly to an A-nucleon system. 

Consider an arbitrary contribution to an irreducible graph involving A nucleons. This 
being only part of the full amplitude, it will in general have C separately connected pieces 
[C = 1 for A = 0,1] C = 1, ...,A — 1 for A > 2.) The straightforward generalization of Eq. 
( |109|) is [0,|ll 

4-A + 2{L-C)+J2VAi. (133) 



v 



Again, because C is bounded from above (C < Cmax), the chiral symmetry constraint 
Aj > still implies a lower bound on the power oi Q, u > Vmin = 4 — A — 2Cmax for strong 
interactions, so that systematic calculations can be carried out. 

With this power counting we can get some insight into few-body forces (those diagrams 
with C = 1). The new forces that appear in systems with more than two nucleons have 
been considered in Refs. [|r^,|T^,|T^,|r^ . The dominant potential, at z/ = 6 — ?>A = Umin, is 
simply the two-nucleon potential of lowest order that we have already encountered in Sect. 
[IV C| . We can easily verify that a three-body potential will arise at z/ = i/mm + 2, a four-body 
potential aX u = fmin + 4, and so on. It is (approximate) chiral symmetry therefore that 
implies that n-nucleon forces VnN are expected to obey a hierarchy of the type 
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{Vi^n+l)N)KVnN) ~ 0((Q/Mqcd)'), (134) 

with (Kiat) denoting the contribution per n-plet. This is a non-trivial consequence of chiral 
symmetry, as there exist non-chiral models that produce large three-body forces. If | (V2Ar) | ~ 
10 MeV as in Sect. [TVTl , we can estimate | (Vg^v) | ~ 0.5 MeV, | (V47V) | ~ 0.02 MeV, and so on. 



This is in accord with detailed few-nucleon calculations using phenomenological potentials. 
It proves instructive to look at the form of the first few terms in the few-nucleon potential 



T^ , p!5|JT7| . We do so with time-ordered perturbation theory, as before. At z/ = Vmin + 2, in 
addition to corrections to the two-nucleon potential, one also finds diagrams involving either 
three nucleons or two pairs of nucleons connected via leading contact interactions and static 
pions. One finds [^ that the various orderings of Fig. [l^(a) add to zero. Diagrams like 
those in Fig. |l^(b,c) give non-vanishing contributions to three- and double-pair potentials. 



However, it is easy to prove ||13| , p!7| that these contributions cancel, to this order in the 
expansion, against contributions from reducible graphs like those in Fig. |T^(d,e) involving 
the leading and the energy-dependent pieces of the two-nucleon potential. This happens 
because a small energy denominator combines with the small energy- dependent OPE to 
produce a result not only of the same order, but exactly opposite to the three- and double- 
pair potentials. This cancellation had been noted before in the case of the TPE 3A^ force 
||130|,|99|, but its model independence and generality are particularly clear in the EFT. Refs. 



107| , |105|| have emphasized that redefining the potential to eliminate energy dependence leads 
to no 3A^ TPE forces of this type at all. Remaining, to this order, are only "true" three-body 
forces generated by the delta isobar, and at z/^m + 3 further terms with similar structure 



arise |TB| , p!7| — see Fig. |TB[ Up to this order, there are no four-nucleon forces, and the 
three-nucleon potential has components with three different ranges: neglecting relativistic 
corrections, 

/ON 

^3 (%■> Qjk) = eiti ■ tk + e2(Jj ■ (Tfctj ■ tk + e-idj ■ {Bi x ak)tj ■ {U x tk) 

— -^^afc ■ qjk [diiti ■ tkffi + tj ■ ikdj) - 2^2*^ ■ {U x tk)ai x dj] ■ Qj 

■K jk 

f2gAY 1 ^ ^ ^ ^ 



k 



X 



ti ■ tkicsQij ■ Qjk + 2ciml) + dtj ■ {ti X tk)aj ■ {qij x q^k) 



+ two cyclic permutations of (ijk), (135) 

with Piipi) the initial (final) momentum of nucleon i, qij = pi — p^ = p- — pj, and wfj = 
qlj +m^. The parameters q and di are given by Eqs. (^) and (|102|) in terms of the Bi from 
Eq. dl), the vrA^A coupling Ha from Eq. (|^, the Di from Eq. ([T00|) , and the NNNA 
contact interaction Dt from Eq. (|99|) . Likewise, the Cj are parameters that have delta and 
shorter-range components Ei, 

ei = El, 62=^2 + ^, 63 = ^3--^. (136) 

9 dm loom 



The TPE part of the potential ( |135|) is determined in terms of vrA^ scattering observ- 



ables. It is not identical to the Tucson- Melbourne potential |^. The latter was built from 
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FIG. 17. Some time-ordered diagrams whose contributions add to nothing in the EFT: (a) TPE 
from the Weinberg- Tomozawa seagull, of which only one ordering is shown; (b,c) "uncorrelated" 
TPE three-nucleon potential and (d,e) iteration of energy-dependent OPE two-nucleon potential. 
Solid lines represent nucleons, dashed lines pions, a heavy dot an interaction in C^^', and a dot 
within two circles an interaction in C^'^'. The same cancellation (b-e) occurs in other sets of 
three-nucleon diagrams corresponding to other orderings. It also takes place in diagrams where 
short-range interactions are substituted for static OPE, and in sets of double-pair diagrams. 
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FIG. 18. Some time-ordered diagrams contributing to the three-nucleon potential V37V in the 
EFT. (Double) solid lines represent nucleons (and/or deltas), dashed lines pions, a heavy dot an 
interaction in C^^', and a dot within a circle an interaction in C^^'. First line corresponds to 



', + 2, second line to z^ 



+ 3, and ". . ." denote v > i'min + ^- All nucleon permutations 



and orderings with at least one pion or delta in intermediate states are included. 



TiN scattering through a particular choice of pion fields for which chiral symmetry is not 
respected term- by-term [|131|| . Although the corresponding on-shell irN amplitude gives the 
same result as Eqs. (^) and (p6D, they differ off shell. The corresponding three-body po- 
tential involves other interactions to guarantee chiral symmetry and to enforce invariance 
under pion-field redefinitions, but unfortunately such terms have not been included in the 
original Tucson-Melbourne potential. After this is done, all commonly used TPE three- 
nucleon potentials agree, except for the numerical values of parameters that depend on the 
way irN data is fitted — see Table [|131|| . The novel one-pion/short-range components of 
the potential (|135|) can be related to pion production in NN collisions. They are starting 
to be explored in conjunction with more conventional treatments of few-nucleon systems 
||138| , |139|| . The purely short-range components can only be determined from few-nucleon 
systems, but to leading order depend only on one undetermined parameter, Dt- (Recall 
that Dt could not be isolated in the z/ = Umin + 3 EFT fit of A^A^ phase shifts.) Relativistic 
corrections neglected above are discussed in Ref. 



65 



TABLE V. Low-energy vrA'^ scattering parameters (with Born terms removed) in a variety of 
existing TPE three- nucleon forces. 



Three-Nucleon Force 


(a — 2m^c) m,r 


bml 


cml 


dml 


Fujita-Miyazawa 132|] 


0.0 
-1.03 
-1.05 

0.0 

-1.87 
-0.51 


-1.15 
-2.62 
-2.29 
-1.20 
-3.82 
-1.82 


0.0 
1.03 
1.05 
0.0 
0.0 
0.0 


-0.29 


Tucson-Melbourne |48,133| 


-0.60 


Brazil 1 134 133 1 


-0.77 


Urbana-Argonne [|135|,136^ 


-0.30 


Texas 0,11 
Ruhr (Pot) [137|| 


-1.12 
-0.48 



VI. PROBES OF FEW-NUCLEON SYSTEMS 



Further constraints on our understanding of few-nucleon systems come from the study 
of processes involving external probes, such as pions, photons, electrons and neutrinos. 
In some cases the calculation in the EFT to a certain order is completely determined by 
parameters already known from other processes; a prediction can be made, accompanied by 
an educated guess about the error resulting from neglecting higher orders. Other cases can 
be used to determine so-far unknown parameters; this is particularly useful in the case of 
neutron and neutral pion parameters that cannot be easily determined in processes with a 
single-nucleon target — see Sect. |ITI|for examples. But even in these cases, a prediction for 
energy dependence is usually possible. 



A. Very low and low energies 

In the EFT with virtual pions integrated out, scattering of pions on few-nucleon systems 
is in itself not very interesting, as pions behave just as heavy particles. More interesting 
are processes involving photons of momenta Q ^ 'R, which we can use as efficient probes of 
momentum distributions. As we have seen, this EFT is equivalent in the two-nucleon sector 
to effective range theory. An extensive literature exists using this technique in the calculation 
of electroweak process involving the deuteron, starting with the pioneering work of Refs. 
l|. 140 . These calculations can be recast in terms of the very- low-momentum EFT in a 



straightforward way. We briefly describe here some of the more modern calculations, carried 
out with explicit but perturbative pions, that is, in the low-momentum EFT. Calculations 
involving the two-nucleon system proceed by separating "irreducible" contributions, which 
are those that do not split in two when cut at any Co vertex. Always contributing is an 
irreducible two-nucleon function, intimately related to the two-nucleon amplitude. Other 
irreducible components include the external probes. Each irreducible part is calculated to the 
same order in the Q/M^j^ expansion, they are combined in the amplitude of the processes 
of interest, and only terms up to the desired power of Q/Mjq^ are kept. These calculations 
typically contain the result of the lower-scale EFT; improvement due to the additional pion 
effects, if any, can be taken as evidence for M^tat > rriT,. 

• Deuteron form factors. A calculation of the deuteron form factors in subleading order 
was carried out in Ref. |p.41 ]. To this order the deuteron is pure S wave, determined by 
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TABLE VI. Results from EFT at leading (LO) and subleading (SLO) order and experimen- 
tal values for deuteron rms charge radius {{r'^)Jf^ ), magnetic dipole moment (fid), and electric 
quadrupole moment (Qe)- 



quantity 



Qe 



LO 



NL+SLO 



expt 



L53 

0.88 




L89 

0.86 (fit) 

0.40 



2.1303(66) 

0.857406(1) 

0.2859(3) 



the ^5*1 parameters Cq + mlC2qm and C2. Interactions with a photon field come, at this 
order, from known one-nucleon terms, from minimal substitution in the C2 term, and from 
one a priori unknown two-body magnetic-interaction counterterm, which can only be fitted 
to data. With Cq + m\C2qm fitted to the deuteron binding and C2 to the phase shifts, the 
results for static moments are in Table |V| ||141[| . The momentum dependence of the charge 
and magnetic form factors can then be predicted. Good convergence and agreement with 
data are attained even at momenta as large as ~ 400 MeV; only the quadrupole moment 
works poorly. Unfortunately no improvement over the effective range expansion is seen. 
Comparison with effective range theory for the radius and EFT with non-perturbative pions 
for the quadrupole moment suggest that the most important of the neglect terms are from 
iterated pion exchange ||141|| . Sub-subleading contributions to the quadrupole moment are 
discussed in Ref . ||142|| . The anapole moment of the deuteron was studied in Ref . ||143|| . 

• np —>■ jd at threshold. To subleading order this reaction proceeds from the ^5*0 state; 
the same set of ^5*0 and ^5*1 parameters Cq + m'^C2qm and C2 appear and are determined by 
phase shifts and deuteron binding. In leading order only one-body currents from magnetic 
coupling appear. At subleading order pion-exchange currents are calculated and found to 
be cutoff dependent; a four-nucleon/one-photon operator appears at the same order and 
absorbs this dependence leaving an undetermined finite part behind, which can then be 
fitted. Results | 144 | are presented in Table |V11| . The parity- violating asymmetry in the 
angular distribution of gamma rays from radiative capture of polarized cold neutrons was 
calculated in Ref. [ |145|| . 

• '-fd —>■ 'jd. Amplitudes with two external photons give rise to interesting predictions. In 
leading non-vanishing order only one-nucleon operators from subleading interactions enter, 
including the Thomson seagull operator and interactions from minimal coupling. At sub- 
leading order pion exchange and C2 appear as well. Scalar and tensor, electric and magnetic 



TABLE VII. Values for the total cross-section a in mb for radiative neutron-proton capture: 
leading (LO) and subleading (SLO) order EFT and experiment (expt). 



LO 



NL+SLO 



expt 



297.2 



334.2 (fit) 



334.2±0.5 
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TABLE VIII. Scalar and tensor electric (ao and 02) and magnetic (/3o and (32) polarizabilities 
of the deuteron in fm^: leading (LO) and subleading (SLO) order EFT 

polarizability LO NL+SLO 

ao (fm^) 0.386 0.595 

Q2 (fm^) -0.062 

/?o (fm3) 0.067 ? 

P2 (fm^) 0.195 ? 



polarizabilities have been computed in Ref. | ]146| | and are shown in Table |V111| . Convergence 



seems consistent with static electromagnetic moments, and results are comparable to the 
effective range expansion. Nucleon polarizability enters as a higher-order contribution, and 
therefore it is unlikely that it can be isolated from the deuteron polarizabilities. Differential 
cross-sections for Compton scattering at photon energies of 49 and 69 MeV were calculated 
in Ref. ||147||. The parameter- free results are in good agreement with existing data, in partic- 



ular at the lower energy. Nucleon polarizabilities from pion loops enter to subleading order 
and are important for this agreement. Tensor polarized scattering was computed in Ref. 



148 1 . At 90° the first non- vanishing contribution comes only from the pion exchange, and 



therefore might allow a sensitive test of the goodness of perturbative pions. 

B. Moderate energies 



The power counting arguments of Sect. [VB| can be easily generalized to the case where 



external particles with momenta Q ~ Mnuc interact with few-nucleon systems. We define 
the kernel K as the sum of irreducible diagrams to which the probes are attached. A generic 
diagram contributing to a full amplitude will consist of irreducible diagrams sewed together 
by states with small energy denominators. Interactions among nucleons occurring before 
or after scattering can be treated as before: iteration of the potential gives rise to the 
wave-function lip) {\ip')) of the initial (final) nuclear state. The full scattering amplitude is 
then 

T ^ {tlj'\K\tlj). (137) 

The power counting ( [L3iJ| ) applies equally well to the kernel K. In practice, it is frequently 
desirable to minimize nuclear wave-function errors by using a high-precision phenomenolog- 
ical potential instead of the only currently available EFT potential |]l6l . That this is a good 
approximation is suggested by a comparison [Q between a simplified version of the EFT 
potential of Ref. [jl6| and the Argonne V18 potential [^, which show agreement in most of 



the wave-function features. 

Before plunging into hard results, let me point out the generic results of this EFT. 



Because of the factor — 2C in Eq. ( p.33| ), we see immediately — in an effect similar to 
few-nucleon forces — that external probes will tend to interact predominantly with a single 
nucleon, simultaneous interactions with more than one nucleon being suppressed by powers 
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of {Q/MqcdY ■ Again, this generic dominance of the impulse approximation is a well-known 
result that arises naturally here. This is of course what allows extraction, to a certain 
accuracy, of one-nucleon parameters from nuclear experiments. A valuable by-product of 
the EFT is to provide a consistent framework for one- and few-nucleon dynamics, where 
few-nucleon process can be used to infer one-nucleon properties. More interesting from the 
purely nuclear-dynamics perspective are, however, those processes where the leading single- 
nucleon contribution vanishes by a particular choice of experimental conditions, for example 
the threshold region. In this case, two-nucleon contributions, especially in the relatively 
large deuteron, can become important. Further examination of the structure of the chiral 
Lagrangian reveals that two-nucleon contributions tend to be dominated by pion exchange. 
Indeed, photons and pions couple to four-nucleon operators only at 0{Q/Mqcd) relative 
to pion-exchange diagrams constructed out of the leading order Lagrangian. Thus power 
counting justifies the "chiral filter hypothesis" that has been put forward to summarize 
some empirical results on electroweak form factors [^. This "pion dominance" ensures 
that two-body contributions from the EFT in lowest orders tend to be similar to those in 
phenomenological models that include pion-exchange currents. 

• 7id —>■ ird. This is perhaps the most direct way to check the consistency of xPT in 
one- and few-nucleon systems. For simplicity, consider the region near threshold. There the 
lowest-order, u = Umm = ~2 contributions to the kernel vanish because the pion is in an S 
wave and the target is isoscalar. The z/ = Urnin + 1 term comes from the (small) isoscalar 
pion-nucleon seagull, related in lowest-order to the pion-nucleon isoscalar amplitude 6^°-'. 
t^ = Vmin + '^i ^min + 'i coutributious couic froui corrcctious to ttN scattering and two-nucleon 
diagrams, which involve besides 6'-°^ also the much larger isovector amplitude fo^^-* . Weinberg 
T3| has estimated these various contributions to the vrd scattering length, finding agreement 



with previous, more phenomenological calculations, which have been used to extract h'^^\ 

The impact of the deuteron process on the determination of h^^' is further discussed in Ref. 

149| |. The size of two-nucleon, one-loop diagrams that appear in next order is estimated to 



be small in the related, double-charge-exchange process |150|. Charge symmetry breaking 



effects are considered in Ref. |p.51|| . 



• np — > 7(i. This offers a chance of a precise postdiction. Here it is the transverse 
nature of the real outgoing photon that is responsible for the vanishing of the lowest-order, 
V = Vmin = —2 contribution to the kernel. The single-nucleon magnetic contributions come 
at u = Vrnin + 1 (tree level), z/ = Prnin + 2 (one loop), etc. The first two-nucleon term is 
one-pion exchange at i/ = Vmin + 2, long discovered to give a smaller but non-negligible 
contribution. There has been a longstanding discrepancy of a few percent between these 
contributions and experiment. At z/ = Vmin + 4 there are further one-pion exchange, two- 



pion exchange, and short-range terms. In Ref. ||152|| the two-pion exchange diagrams in 
the deltaless theory were calculated and resonance saturation used to estimate the other 
i^ = ^min + 4 terms. With wave-functions from the Argonne V18 potential |Q and a cut-off 



A = 1000 MeV, the excellent agreement with experiment shown in Table PC| was found ||152 
The total cross-section changes by only 0.3% if the cut-off is decreased to 500 MeV. 

• 7(i ^ 7(i. External photons couple to the kernel only in subleading order, v = Umm + 
1 = —1, via the one-body Thomson seagull. In next order, z/ = Umin + 2, more interesting 
effects are present, such as a pion loop which contributes to the nucleon polarizability. There 
are also two-nucleon contributions from pion exchange, which can and need to be computed, 

69 



TABLE IX. Values for various contributions to the total cross-section a in mb for radiative 
neutron-proton capture: impulse approximation to z^ = Vmin + 4 (imp), impulse plus two-nucleon 
diagrams at v = Umin + 2 (imp+tnO), impulse plus two-nucleon diagrams up to z^ = Vmin + 4 
(imp+tn), and experiment (expt). 

imp imp+tnO imp+tn expt 

305.6 321.7 336.0 334+05 



if one seeks to extract information about the neutron polarizabiUty. Work is in progress ||153| 
to calculate the differential cross-section at 49, 69, and 95 MeV. Preliminary results are in 
good agreement with the existing Illinois data at the lowest two energies, and a prediction 
will be made for the highest energy, currently being measured and analyzed in Saskatoon. 

• 7(i — > Tx^d. This reaction offers the possibility to test a prediction arising from a 
combination of two-nucleon contributions and the neutral-pion single-neutron amplitude. 
The differential cross-section for a photon of momentum k to produce a pion of momentum 
q is, at threshold, [{k / q){da / dVL)]q=Q = 8\Ed\'^/3. Ed has been studied up to z/ = z/^m + 3 
with the delta integrated out in Ref. [15^]. There are two classes of contributions, according 



to whether the external light particles interact with one or with both nucleons. The one- 
nucleon part of the kernel is given by the same z/ = Umin, ^min + ^, ■■■ mechanisms described in 
Sect, m, with due account of P waves and Fermi motion inside the deuteron. The neutrality 
of the outgoing 5'-wave pion ensures that the leading u = —2 = z/^m terms vanish. The first 
two-nucleon part of the kernel appears at z/ = z/^m + 2; it comes from a virtual charged pion 
photoproduced on one nucleon which rescatters on the other nucleon with charge exchange. 
These contributions are actually numerically larger than indicated by power counting due 
to relatively large deuteron size. Smaller two-nucleon terms appear at z/ = Umin + 3 from 



corrections in either nucleon. Results for E^ up to z/ = Umin + 3 ||154|| are shown in Table 



^ They correspond to the Argonne V18 potential ||8^ and a cutoff A = 1000 MeV. Other 



realistic potentials and cutoffs from 650 to 1500 MeV give the same result within 5%. The 



chiral potential of Sect. |IVC| is more cumbersome to use, but it has been verified that it gives 
J^ = J^min + 2 results that are similar to other realistic potentials. We see that two-nucleon 
contributions seem to be converging, although more convincing evidence would come from 



next order, where loops appear. A model-dependent estimate [|155|| of some u = Vmin + 4 
terms suggests a 10% or larger error from neglected higher orders in the kernel itself. The 
single-scattering amplitude depends on amplitude for 772 -^ n^n, Eo^^n'^n), in such a way 
that Ed ~ -1.79-0.38(2.13-Eo+(7r°n))inunitsoflO-V^7r+- Thus, sensitivity to Eo+(7r°n) 
survives the large two-nucleon contribution at z/ = Umin + 2. 

We see that working within the effective theory yields a testable prediction, E^ = 
— (1.8 ± 0.2) • 10^^/171.^+. It is remarkable that for this process xPT gives results that are 
somewhat different from tree-level models of the type traditionally used in nuclear physics. 



For example, the models in Ref. ||156|| predict the threshold cross-section about twice as large 



as xPT. Most of the difference comes from one-nucleon loop diagrams: tree-level models tend 
to differ from xPT mostly by having a smaller Eo^{7i^n), which increases \Eii\. A test of 
this prediction is thus an important check of our understanding of the role of QCD at low 
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FIG. 19. Reduced cross-section o"^ = {k/q)a in /xb for neutral pion photoproduction as 
function of the photon energy in MeV. Threshold is marked by a dotted hne. Squares are data 
points from Ref . |157| and the star at threshold is the xPT prediction of Ref . [ 154 1 . Figure courtesy 
of U.-G. Meifiner. 



energies. Such a test was recently carried out at Saskatoon [ p.57|| . The experimental results 
for the pion photoproduction cross-section near threshold are shown in Fig. |19|, together 
with our xPT prediction at threshold. Inelastic contributions have been estimated in Ref. 
1 157 1 and are smaller than 10% throughout the range of energies shown. At threshold, Ref. 
| 157|| finds E^ = — (1.45±0.09)-10~^/m,r+- While agreement with xPT to order z/ = i'min + ^ 
is not better than a reasonable estimate of higher-order terms, it is clearly superior to tree- 
level models. This is compelling evidence of chiral loops. Further test will come from an 
electroproduction experiment currently under analysis in Mainz |[158| . 

• Axial currents. The long-range two-nucleon contributions to the axial current come 
from one-pion exchange in first non-vanishing order and from two-pion exchange as a first 
correction. They have been evaluated in Ref. 
pp — ^ de^Pe is found |^ in agreement with previous calculations. 



14 1 . The related proton-burning process 



pp 



ppn dose to threshold. This reaction has attracted a lot of interest because 



TABLE X. Values for E^ in units of W~^/m^i 
V = ^min + 3, two-nucleon kernel (2A^) aX u = v„ 



from one-nucleon contributions (lA^) up to 
;„ + 2 and at z^ = Vmin + 3, and their sum 



(liV + 2iV). 



liV 




2A^ 




1N + 2N 


V < J^min + 3 


l" = J^min + 2 




I^ = I^min + 3 


V < fmin + 3 


0.36 


-1.90 




-0.25 


-1.79 
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FIG. 20. Cross-section for pp — > ppn^ in fih/sv as function of the pion momentum rj in units of 
mvr for two NN potentials (Argonne V18 and Reid93) and two parametrizations of the isoscalar 
pion-nucleon amphtude (ste and cl). 



of the failure of standard phenomenological mechanisms in reproducing the small cross- 
section observed near threshold. It involves larger momenta of 0{y/rftPmN), so the relevant 
expansion parameter here is the not so small {171^^/171^)2 . This process is therefore not a 
good testing ground for the above ideas. But (m^/mTv)^ is still < 1, so at least in some 
formal sense we can perform a low-energy expansion. In Ref. |^^ the chiral expansion was 
adapted to this reaction and the first few contributions estimated. (Note that — contrary to 
what is stated in Ref. ||159|| — momenta ~ ^JrrQmjq do not necessarily imply a breakdown 
of the non-relativistic expansion, a.s p'^/m% ~ (m^/miv)(p^/'n^Af)-) Again, the lowest order 
terms all vanish. The formally leading non-vanishing terms — an impulse term and a similar 
diagram from the delta isobar — are anomalously small and partly cancel. The bulk of the 
cross-section must then arise from contributions that are relatively unimportant in other 
processes. One is isoscalar pion rescattering for which two sets of xPT parameters were 
used: "ste" from a sub-threshold expansion of the vrA^ amplitude [^] and "cl" from an one- 
loop analysis of threshold parameters ^^ . Others are two-pion exchange and short-range 
ttNNNN terms, which were modeled by heavier-meson exchange: pair diagrams with a 
and u exchange, and a irpu coupling, among other, smaller terms 
Argonne V18 [H and Reid93 [[M|] - 
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Two potentials — 
are shown in Fig. ^ together 



were used. Results 

with lUCF ||162|| and Uppsala [|163|| data. Other xPT studies of this reaction can be found 
in Ref. [|164|| . The situation here is clearly unsatisfactory, and presents therefore a unique 
window into the nuclear dynamics. Work is in progress, for example, on a similar analysis 
for the other, not so suppressed channels — > dn^, -^ prni^ ||165| . 
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VII. MANY-NUCLEON SYSTEMS 

The density of a many-body system sets the scale of typical particle momenta Q. At 
sufficiently low density such that Q < N, we have seen that the EFT can be formulated in 
terms of contact interactions only, and should be perturbative. There is a vast literature on 
dilute systems (see, e.g. Ref. |2^). Most of its results can be rederived in an EFT language. 



For example, Ref. [16C] has reexamined some of the higher corrections to the energy density 



of a dilute boson gas; and Ref. ||167|| has studied pairing in a dilute Fermi system, producing 
analytical expressions relating the pairing gap, the density, and the energy density to the 
scattering length. Only exploratory work has so far been carried out regarding the more 
challenging, higher densities where Q/N ^ 1. This should not be surprising in view of 
the non-trivial aspects that even the three-body system bears, as discussed in Sect. |V A. 



Attempts to look at nuclear matter at equilibrium density and at finite nuclei have been 
phenomenological; within its limited scope, this pioneering work has nevertheless uncovered 
remarkable systematics. 

It was very quickly realized that chiral Lagrangians could be related to Walecka-type 
models. Under the assumptions that pion and spin effects average to very little in heavy 
nuclei, the main interactions should be given by contact interactions. Under the further 
approximation of mean field, two four-nucleon contact interactions can indeed be related 



to the original Walecka model, with parameters of the correct order of magnitude ||168 
More elaborate versions of the Walecka model can be shown to require further contact 
interactions. The more extensive analysis of Ref. |p,69|| has shown a remarkable regularity 



in the size of parameters in a relativistic point-coupling model that includes several four-, 
six-, and eight-nucleon interactions, is treated in Dirac-Hartree mean-field, and is fitted to a 
large number of nuclear data. This regularity was confirmed and extended to more complete 



versions of the same model |p.7CI|| , and was shown to also hold in non-relativistic, Skyrme-force 
models [|171|| . It was found that a 2n-nucleon field operator without derivatives contributes 
~ {4:Po/F^Mqcd)^~^Mqcd to the energy per particle at equilibrium density po- This is 
tantalizingly reminiscent of our previous naturalness assumptions. It has also been noted 
| 171| | that the cancellation between large scalar and vector densities that is characteristic of 



relativistic point-coupling models is incorporated automatically in a smaller contribution in 
the non-relativistic models. Although these models do not include all operators consistent 
with symmetries to a certain order in a definite power counting, additions of new interactions 
with natural parameters could change the precise values of all existing parameters but not 
their orders of magnitude. This regularity is a nontrivial result because at mean-field level 
effective interactions absorb long distance many-body effects from ladder and ring diagrams 
1 1 72| | . It implies a convergent density expansion for mean-field contributions to nuclear 



matter, as the expansion parameter is Apo/F^MqcD ~ 0.2 [|T^. For an attempt at an EFT 



expansion around po, see Ref. [ [I73|| . This all hints that a controlled EFT expansion for finite 
density nuclear systems is possible. 
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VIII. OUTLOOK 

We have gone a long way in identifying the essential ingredients of the EFT appropriate 
for nuclear physics. 

• In the two-nucleon system, the way fine-tuning surfaces in contact interactions, the 
role of a potential and of a dimeron field, and the virtues and limitations of various 
regulators have been understood. As a consequence, systematic calculations can and 
have been performed, in some cases producing simple analytical results. Two-nucleon 
observables work well at low energies, but not yet at the level of phenomenological 
potential models. 

• In the three-nucleon system, we have learned about the subtlety of renormalization in a 
non-perturbative context and started exploring the rich features of phenomenology and 
three-body forces. Universal features — such as the Phillips line — arise as consequence 
of the renormalization group. 

• Nuclear processes involving external particles have been treated consistently with one- 
nucleon xPT, providing a way to access neutron properties. In most cases results 
consistent with conventional models have been obtained, but in one case there was 
a j^rediction in disagreement with these models. This has stimulated experimental 
efforts which unequivocally favor the EFT. 

• Evidence for naturalness has been unearthed in many-nucleon systems. 

When these developments are taken together, a simple physics picture of nuclei emerges 
from QCD. Low-momentum nucleons can be thought of as point particles surrounded by (i) 
an inner cloud which is dense but of short range 0{1/Mqcd), and (ii) an outer cloud made 
out of pions which has long range 0(1/M„„c) but is sparse. If a few non-relativistic nucleons 
are put together, each nucleon is unable to distinguish details of the others' inner clouds. 
Interactions in the inner region can be expanded in delta functions and their derivatives, 
of progressively smaller importance. The outer cloud yields interactions which are non- 
analytic, but chiral symmetry guarantees that the chance of finding n pions in the air at the 
same time decreases rapidly with n. The systematic theory built on this simple picture has 
a finite number of interactions at any desired accuracy. It in fact explains many regularities 
that have been discovered in previous phenomenological treatments but whose origins are 
otherwise mysterious. These include the decreasing importance of many-pion exchanges, 
many-nucleon forces, and higher isospin classes. 

This picture is not in contradiction with potential-model experience. In fact, the method 
of EFT has already provided quantitative input for more standard analyses that could 
not have been obtained before by other means. The chiral TPE two-nucleon force and 
simultaneous pion-gamma exchange have been successfully incorporated in the Nijmegen 
phase shift analysis; chiral TPE three-nucleon forces have been used to correct existing 
models; and, at the same time, the effects of the shorter-range three-nucleon forces on 
discrepancies remaining in conventional approaches are starting to be explored. From a 
conservative point-of-view, the EFT method can be considered a refinement of other hadronic 
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approaches in that it brings to an otherwise free-fantasy context the constraints of QCD 
symmetries. 

There is no denying that much is still to be done in raising this pre-teen EFT formulation 
of nuclear physics into a fully mature subject. One major source of malnourishment has been 
the so-far elusive determination of the scale where pion effects cease to be perturbative, and 
the scope of their resummations. Yet, this is a problem that can and will be resolved. A 
coherent description of few- and many-nucleon systems is still green but we can already 
glimpse the shape it is going to take. 
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